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ABSTRACT 

In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane 
telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the ob¬ 
servations, reductions and measurements, and data products from the spectroscopic campaign. The reduced 
spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum 
and emission-line components. We present light curves of broad emission lines and the AGN continuum, and 
measurements of the broad H/3 line widths in mean and root-mean square (rms) spectra. For the most highly 
variable AGNs we also measured broad H i line widths and velocity centroids from the nightly spectra. In 
four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H/3 width 
and luminosity, demonstrating that the broad-line region “breathes” on short timescales of days to weeks in re¬ 
sponse to continuum variations. We also find that broad H/3 velocity centroids can undergo substantial changes 
in response to continuum variations; in NGC 4593 the broad H/3 velocity shifted by ~ 250 km s _I over a 
one-month duration. This reverberation-induced velocity shift effect is likely to contribute a significant source 
of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary 
orbital motion. We also present results from simulations that examine biases that can occur in measurement 
of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting 
noise. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert— techniques: spectroscopic 


1 Department of Physics and Astronomy, 4129 Frederick Reines Hall, 
University of California, Irvine, CA, 92697-4575, USA; barth@uci.edu 

2 Physics Department, California Polytechnic State University, San Luis 
Obispo, CA 93407, USA 

3 Department of Physics and Astronomy, University of California, 
Riverside, CA 92521, USA 

4 Department of Astronomy, University of California, Berkeley, CA 
94720-3411, USA 

5 Lick Observatory, P.O. Box 85, Mount Hamilton, CA 95140, USA 

6 Department of Astrophysical Sciences, Princeton University, Prince¬ 
ton, NJ 08544, USA 

7 Department of Physics and Astronomy, University of California, Los 
Angeles, CA 90095-1547, USA 

8 Department of Physics, University of California, Santa Barbara, CA 
93106, USA 

9 Texas Tech University, Physics Department, Box 41051, Lubbock, TX 
79409-1051, USA 

10 Jet Propulsion Laboratory, California Institute of Technology, 4800 
Oak Grove Boulevard, Pasadena, CA 91109, USA 

11 Astronomy Program, Department of Physics and Astronomy, Seoul 
National University, Seoul 151-742, Republic of Korea 

12 Nucleo de Astronomia de la Facultad de Ingenierfa, Universidad 
Diego Portales, Av. Ejercito Libertador 441, Santiago, Chile 

13 Department of Astronomy and Center for Galaxy Evolution Research, 
Yonsei University, Seoul 120-749, Republic of Korea 

14 Department of Statistics, The University of Auckland, Private Bag 
92019, Auckland 1142, New Zealand 

15 Astrophysics Science Division, NASA Goddard Space Flight Center, 
MC 661, Greenbelt, MD 20771, USA 

16 Joint Space-Science Institute, University of Maryland, College Park, 
MD 20742, USA 

17 Department of Astronomy, University of Wisconsin, Madison, WI 
53706, USA 

18 FONDECYT Postdoctoral Fellow, Departamento de Astronomia, 
Universidad de Concepcion, Chile 


19 School of Physics & Astronomy, University of Southampton, 
Southampton SOI7 1BJ, UK 

20 Department of Physics and Astronomy, Aarhus University, Ny 
Munkegade, DK-8000 Aarhus C, Denmark 

21 Dept, of Physics and Physical Science, University of Nebraska at 
Kearney, Kearney, NE 68849, USA 

22 Center for Cosmology and Astro-Particle Physics, The Ohio State 
University, Columbus OH 43210, USA 

23 Department of Astronomy, The University of Texas, Austin, TX 
78712-0259, USA 

24 Astronomy Department, Yale University, New Haven, CT 06510, 
USA 

25 George P. and Cynthia Woods Mitchell Institute for Fundamental 
Physics and Astronomy, Department of Physics and Astronomy, Texas 
A&M University, College Station, TX 77843-4242, USA 

26 Deceased 12 December 2011 

27 Hubble Fellow 


2 


1. INTRODUCTION 

Broad emission lines are a hallmark property of active 
gal actic nuclei (AGNs), first highlighted in pioneering work 
by ISevferti (119431) . Originating from the broad-line region 
(BLR), they serve as important diagnostics of physical condi¬ 
tions within the AGN as well as dynamical tracers providing 
information about the mass of the supermassive black hole at 
the core of the AGN central engine. Unfortunately, even for 
the nearest AGNs the BLR is too small to be directly resolved 
by any current observational capability. For example, in the 
well-studied nearby Seyfert 1 galaxy NGC 5548, the radius 
of the H.3-emitting portion of the BLR is ~ 4-20 light-days 
dBentz et al.ll2007t) . corresponding to an angular radius of just 
9-45 microarcseconds. 

Although the spatial structure of the BLR cannot be directly 
resolved, the time-variable nature of AGN emission makes it 
possible to resolve the BLR’s size and structure in the time do¬ 
ma in via reverberation mapping ( iBlandford & McKeell 19821 
lPetersonll2001l) . Ionizing photons from the AGN’s central en¬ 
gine are reprocessed by BLR gas into emission-line photons. 
From the vantage point of a distant observer, random varia¬ 
tions in the flux of the AGN continuum will be followed by 
corresponding variations in the fluxes of emission lines from 
the BLR, with a time delay that depends on the light-travel 
time across the BLR. Cross-correlation of emission-line light 
curves against light curves of the AGN continuum gives a 
measure of the BLR radius iblr, under the assumption that 
the size of the continuum emission region is much smaller 
than jblr. Reverberation mapping thus requires measurement 
of light curves of the variable AGN continuum and emission 
lines, with a cadence sufficient to resolve the temporal fluc¬ 
tuations adequately. When /blr is combined with a veloc¬ 
ity parameter measured from the width of the broad line, a 
so-called virial mass estimate can be derived for the central 
black hole. Virial masses have been determined through re¬ 
verberation mapping for several dozen low-redshift AGNs, 
and these measurements form the foundation of secondary 
“single-epoch” virial mass methods, which provide nearly all 
of the empirical information currently available on the cosmo¬ 
logical growth history and redshift evolution of supermassive 
black holes (for a comprehensive review, see lShenll201 3l) . 

A major goal for reverberation mapping observations is the 
detection of velocity-resolved emission-line variations — that 
is, variations in the shapes of broad emission-line profiles, 
not merely variations in total emission-line flux. Line profile 
variability conveys information about the kinematics of BLR 
gas. Measurements of profile variability can provide a means 
to discriminate between rotationally dominated, inflowing , 
or outflowing BLR kinematics (e.g., I Welsh & HornellT99l]) . 
while high-cadence and high signal-to-noise ratio (S/N) ob¬ 
servations can potentially reveal_a wealth of detail about BLR 
structure and kinematics (iHorne et al.1120041) . 

It is a nontrivial challenge to obtain observations with suf¬ 
ficiently high cadence and high S/N suitable for measurement 
of broad-line profile variations. Consequently, many of the 
largest reverberation-mapping programs have focused largely 
on measurement of fl ux variations and derivation of black hole 
virial masses (e.g., iKasni et ahl 120001 iPeterson et al.l 120041 
and references therein). Over the past decade or so, various 
observing programs have succeeded in measuring velocity- 
resolved variability signatures in broad H/3; an early exam¬ 
ple was a sparsely sampled program at t he Hobby-Eberly 
telescope targeting Mrk 110 (IKollatschnyl 120031) . More re¬ 


cently, densely sampled observing campaigns have led to 
rapid progress. Using data from the 2008 Lick AGN Mon¬ 
itoring Project campaig n, a 64-night program carr ied out at 
the Lick 3 m telescope, iBentz et ahl (2 0081 2009b ) detected 
a strong gradient in the reverberation lag across the H/3 line 
profile in Arp 151 (Mrk 40), and the campaign provided some 
degree of velocity resolution of the H/3 reverberation lag in 
five additional AGNs. Similarly, intensive monitoring cam¬ 
paigns carried out at the MDM Observatory have recently 
yield ed impressive detections of velocity-res olved reverbera¬ 
tion ( Denney et a l.!'2009b; i Grier et a l.li 2013 a). The collective 
results of these programs suggest a surprising diversity in the 
kinematic properties of AGNs, with different objects exhibit¬ 
ing behavior consistent with rotation, infall, or outflow as the 
dominant kinematic behavior in the BLR. 

These high-cadence monitoring programs at Lick and 
MDM have spurred the development of new methods for anal¬ 
ysis of reverberation data. The traditional approach of deriv¬ 
ing black hole masses by applying a simple recipe to derive 
the “virial product” [VP = /-blr(AV) 2 /G, where /blr is de¬ 
termined from the reverberation lag and AV is a measure of 
the broad-line width] suffers from several shortcomings. In 
particular, the “BLR radius” measured from the reverberation 
lag gives a highly oversimplified measure of BLR structure, 
since the BLR actually spans a very broad range in radius. 
Similarly, a single measure of BLR velocity derived from an 
emission-line width is a gross oversimplification of the com¬ 
plex velocity structure of a real BLR. More worrying is the 
possibility that the BLR dynamics might be affected by out¬ 
flows or other nongravitational motions, in which ca se sim¬ 
ple vi r ial mass estimates mig ht be entirely misleading dKrolikl 
120011) . iPancoast et al.l (1201 ll) proposed a different approach to 
reverberation data analysis, in which black hole masses can be 
determined via a forward-modeling procedure that computes 
the response to continuum variations of a physically moti¬ 
vated dynamical model for BLR clouds. By exploring a broad 
parameter space of BLR dynamical models and comparing the 
model predictions directly with the observed time series of 
broad emission-line profiles, direct constraints can be placed 
on the black hole mass, BLR inclination and opening angle, 
radial emissivity profile, and other parameters. This approach 
has been a pplied to several AGNs from our high-cadence Lic k 
programs dBrewer et alJl20lH IPancoast et all I2012L12014bl) . 
iLi et all (120131) presented applications of an independent im¬ 
plementation of a similar method to archival data from several 
observing campaigns. 

Although reverberation mapping has seen dramatic 
progress over the past several years stemming from these 
high-cadence monitoring campaigns, there are still only a 
handful of AGNs having time-series data with sufficient sam¬ 
pling and S/N for these dynamical modeling methods to 
be optimally successful. In order to obtain more velocity- 
resolved reverberation data suitable for dynamical modeling, 
we embarked on a new reverberation mapping program in the 
spring of 2011. Spectroscopic observations were conducted 
at the Lick 3 m Shane telescope, and U-band photometric 
monitoring was carried out using several queue-scheduled and 
robotic telescopes. This paper describes the spectroscopic 
campaign and presents measurements from the Lick data. 

A new aspect of our data analysis for this 2011 program 
is the use of spectral fitting and decomposition techniques to 
isolate different portions of the variable spectra. In the past, 
nearly all emission-line light curves for reverberation map- 
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ping were measured by subtracting a local, linear continuum 
model underlying strong emission lines such as H/3 in or¬ 
der to isolate the emission-line flux. This method is gener¬ 
ally effective for measuring integrated light curves of strong 
emission lines having high variability amplitude. However, 
simple linear continuum subtraction is unable to remove host- 
galaxy starlight features that can affect the broad-line profiles 
and widths, and it cannot disentangle blended emission lines. 
The alternative approach of fitting models to deblend com¬ 
ponents such as Fe II from H (3 had seldom been applied to 
nightly spectra from reverberation campaigns (e.g. lBian et al.l 
l201d ). Our initial experiments with spectral decomposition 
using observations of Mrk 50 from this program established 
that multicomponent fits made it possible to measure light 
curves for He II A4686 even when the line was heavily diluted 
by starlight, and to detect variabili ty and measure reverbera - 
tion lags of the optical Fe II ble nds dB arth et a l.ll201 1 bl f20 1 3h . 
Using a similar fitting method. iPark et al.l (1201 2 ah carried out 
decompositions of spe ctra from the 2008 Lick AGN Monitor¬ 
ing Project campaign ( Bentz et~akll2009bl) . This work demon¬ 
strated clear improvements in determination of the mean and 
root-mean square (rms) line profiles for H i after removal of 
the host-galaxy starlight and other blended line components 
from the spectra. In this paper, we present a complete de¬ 
scription of our fitting method and results from applying the 
spectroscopic decomposition procedure to our full 2011 mon¬ 
itoring sample. 

The outline of this paper is as follows. Section |2]describes 
the sample selection, and §0 describes the spectroscopic ob¬ 
servations carried out at Lick. Details of the data reduction, 
calibrations, and the fitting procedure applied to the blue-side 
spectra are given in $4] Sections [5]and|6]present the measure¬ 
ments of emission-line light curves, broad H/3 line widths, 
and velocity centroid variations. Section [7] examines the an¬ 
ticorrelation between broad H/3 width and luminosity in the 
single-epoch spectra of highly variable sources. In $8]we de¬ 
scribe how broad-line velocity shifts induced by asymmetric 
emission-line reverberation can potentially mimic the obser¬ 
vational signature of orbital motion in a binary black hole sys¬ 
tem, and we consider the implications of broad-line reverber¬ 
ation for spectroscopic binary black hole searches. Section [9] 
presents a summary of our results and conclusions. In Ap¬ 
pendix [A] we display spectra of additional AGNs that were 
observed during the first few nights of our program as we fi¬ 
nalized our target list for monitoring. Appendices [B] and [Cl 
present simulations designed to clarify the relationship be¬ 
tween line widths measured from mean and rms spectra, and 
to demonstrate potential sources of bias that can afflict broad¬ 
line widths measured from rms spectra. 

2. SAMPLE SELECTION 

The major science goal for this program was to mea¬ 
sure velocity-resolved reverberation signals in H/3 for bright 
Seyfert 1 galaxies, and we primarily selected targets having 
strong broad H/3 emission from the Sloan Digital Sky Survey 
(SDSS) Data Release 7 archive (lAbazaiian et al.l l20Q9t) and 
other AGN catalogs in the literature. Our reverberation cam¬ 
paign began in late March and continued through the mid¬ 
dle of June. This set a requirement for the sample to con¬ 
tain AGNs distributed over a broad range of right ascension 
such that the entire sample could be observed over nearly 
the full duration of the program. Targets were selected to 
have apparent magnitudes of V <17 so that high-S/N spec¬ 
tra could be obtained in short exposures (10-40 min). We 


chose AGNs for which estimated H/3 reverberation lags were 
< 25 days based on either the radius-luminosity relation from 
iBentz et ahi (l2009al) or previous reverberation measurements. 
Finally, the standard configuration of the Kast spectrograph 
employs a 600-line grism in the blue camera and a dichroic 
with cutoff wavelength 5500 A, and the blue-camera through¬ 
put drops steeply at A > 5400 A. This imposes a practical 
redshift limit of z < 0.08 in order for [O III] A5007 to fall on 
the high-sensitivity portion of the blue camera. This line is 
used as the reference for flux calibration of the spectral time 
series. 

Two targets, Mrk 817 and Zw 229-015, were included in the 
sample in order to coordinate with a Spitzer infrared monitor¬ 
ing program taking place in 2011 (Spitzer program ID 70119; 
PI V. Gorjian). Zw 2 29-015 was also a Kepler monitoring 
target during this time dMushotzkv et alJl201 ltlCarini & Rvld 

l2012t lEdelson et al.ll20l4. 


We assembled an initial set of 15 high-priority targets plus 
several alternates which included narrow-line Seyfert 1 galax¬ 
ies, AGNs that were somewhat fainter than our primary tar¬ 
gets, and candidate AGNs for which no high-quality recent 
spectra were available in the literature. During the first few 
usable nights of the campaign, we found that three of our 
top-priority targets selected from SDSS had transformed into 
Type 1.9 Seyferts with essentially no broad H/3 emission re¬ 
maining (see Appendix lAl for details). Several alternate tar¬ 
gets were then observed and we chose the most promising 
ones to fill the gaps in the sample’s right ascension range. Our 
final monitoring sample is listed in Table [Q and AppendixlAl 
presents spectra of the alternate and discarded targets. 

Some of these AGNs have been observed in pre¬ 
vious reverberati o n programs, including Mr k 40 (Arp 


151; IBentz et ail 12008 


1990: Santos-Lleo et al 


1998 jDennevetal 


I. 2010) 

Denne 


I2009bl) . Mrk 279 dMaoz et al 
120011) . Mrk 817 dPeterson et a~ 


NGC 4593 (iDietrich et all 11994 


Onken_etalj 2003 iDennev et al.l 12006 ). and Zw 229-015 
dBarth et al.ll201 la ). Revised lag measurements of Mrk 279, 
Mrk 817, and NGC 4 593 were al so given in the com pilation of 
iPeterson et al.l (120041) . Recentlv. lWang et~aXl (120141) presented 
results from a reverberation campaign that included Mrk 486 
and Mrk 493 as well. New observations of previously well- 
observed targets such as Mrk 40 can potentially test whether 
velocity-resolved reverberation lags remain constant or vary 
over durations comparable to the dynamical timescale of the 
BLR. 


Owing to the combination of criteria given above, the ob¬ 
served sample is not complete or unbiased in terms of black 
hole mass, luminosity, broad-line width, or any other funda¬ 
mental or observed property. This is a recurring issue for all 
AGN reverberation mapping campaigns to date, and should be 
kept in mind when extrapolating reverberation results to de¬ 
rive single-epoch masses for AGNs at high luminosities or Ed¬ 
dington ratios not probed by the reverberation- mapped sam- 
ple. This point h as been emphasized bv lRichards et al.l (1201 ll) 
and IShenl (120131) in the context of examining biases inherent 
in single-epoch masses derived from the C IV line. 


3. OBSERVATIONS 


This project was allocated 69 nights at the Lick 3 m Shane 
telescope, distributed between 2011 March 27 and June 13 (all 
dates are UT). Interspersed with our nights were occasional 
gaps, mostly during or close to full moon, for other projects 
such as exoplanet search programs. 

Observations were conducted using the Kast double spec- 
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TABLE 1 

Sample Properties and Observation Parameters 


Galaxy 

Alt. Name 

Z 

D l 

(Mpc) 

Ay 

(mag) 

Slit PA 
(deg) 

Aobs 

Mean Sampling 
(days) 

Median 

Airmass 

Median 

S/N 

Mrk 40 

Arp 151 

0.0211 

92.3 

0.039 

90 

39 

1.97 

1.07 

72 

Mrk 50 


0.0234 

102.5 

0.044 

180 

55 a 

1.8 l c 

1.29 

77 

Mrk 141 


0.0417 

185.2 

0.028 

90 

36 

1.79 

1.16 

107 

Mrk 279 

PG 1351+695 

0.0305 

134.3 

0.044 

120 

34 

2.30 

1.19 

106 

Mrk 486 

PG 1535+547 

0.0389 

172.4 

0.040 

90 

27 

2.28 

1.08 

94 

Mrk 493 


0.0313 

137.9 

0.068 

70 

32 

2.39 

1.08 

88 

Mrk 504 

PG 1659+294 

0.0359 

158.7 

0.135 

60 

36 

2.18 

1.05 

89 

Mrk 704 


0.0292 

128.5 

0.079 

45 

38 

1.68 

1.15 

70 

Mrk 817 

PG 1434+590 

0.0315 

138.8 

0.019 

90 

27 

2.42 

1.08 

84 

Mrk 841 

PG 1501+106 

0.0364 

161.0 

0.082 

45 

35 

1.82 

1.14 

94 

Mrk 1392 


0.0361 

160.1 

0.125 

45 

38 

2.02 

1.24 

98 

Mrk 1511 

NGC 5940 

0.0339 

149.7 

0.112 

45 

40 

2.00 

1.21 

72 

NGC 4593 

Mrk 1330 

0.0090 

39.0 

0.068 

45 

43 

1.85 

1.44 

121 

PG 1310-108 

II SZ 10 

0.0343 

151.5 

0.143 

180 

35 

2.30 

1.58 

50 

Zw 229-015 


0.0279 

122.6 

0.198 

56.5 

29 b 

2.72 c 

1.07 

49 


Note. — Redshifts are taken from NED. Luminosity distances ( Dl ) are calculated from redshifts assuming a 
WMAP9 cosmology with Ho = 69.7 km s -1 Mpc -1 , Q-m = 0.281, and Oa = 0.7185 <Hinshaw et alJ[20l3h . using the 
IWrighll <20061) calculator. Galactic extinctions (Ay) are from NED and are based on ISchlaflv & Finkbeineil <201 ID . as¬ 
suming Ry = 3.1. N 0 b s is the total number of spectroscopic observations of each AGN. Mean sampling is the mean time 
interval between successive observations, in days. (The median sampling interval is 1.0 days for all AGNs.) The median 
S/N gives the median S/N per pixel between rest wavelengths 4600 and 4700 A in the extracted blue-side spectra, for 
the full series of observations of each AGN. 

a Twelve of the 55 spectra of Mrk 50 were obtained before the start of our main campaign during nights assigned to 
other projects. 

b Three of the 29 spectra of Zw 229-015 were obtained after the end of our main campaign during nights assigned to 
other projects. 

c Mean sampling for Mrk 50 and Zw 229-015 was calculated based on observations taken during the main campaign. 


trograph dMiller & Stonelll99 3). On the blue side (spatial 
scale 0."43 pixer 1 ), we used a 600 lines mm -1 grism cov¬ 
ering 3440-5515 A at 1.02 A pixel -1 . On the red side (spatial 
scale 0"78 pixel -1 ), we used a 600 lines mm -1 grating blazed 
at 7500 A, giving coverage of 5410-8200 A at 2.35 A pixel -1 . 
The d55 dichroic was e mp loyed to split the light between 
the blue and red cameras^ A slit width of 4" was adopted 
in order to mitigate slit losses while also preserving moder¬ 
ate spectral resolution. Calibration frames taken each after¬ 
noon included bias exposures, dome flats, and wavelength- 
calibration lamps including He, Ne, Ar, Hg, and Cd. The red- 
side CCD exhibits a steep change in bias level at the short- 
wavelength end, while the blue-side detector has extremely 
clean and flat bias structure with no discernible persistent fea¬ 
tures. 

The red-side spectra are affected by strong fringing at A > 
7000 A, but it can be removed fairly well by flat-field divi¬ 
sion if a dome-flat exposure is taken at the same telescope 
position and spectrograph rotation angle as the object. While 
the dome-flat exposures are very short (1-2 s), this requires 
additional overhead time for the dome to be rotated or even 
capped (if the target is close to the zenith). In order to maxi¬ 
mize the on-source exposure time during the nights, we opted 
to use only afternoon dome-flat exposures rather than observ¬ 
ing dome flats at the position of each AGN during the nights. 
This trade-off compromises the data quality of the red-side 
exposures at wavelengths longer than 7000 A, but it does not 
affect the blue-side data which are the main focus of our re¬ 
verberation measurements. 

Exposure times for each AGN were between 300 and 2400 

- 8 In our 2008 Lick campaign IBentz et~aOt2009bl) , we used only the red 
arm of the Kast spectrograph because the blue-side CCD had recently failed 
and was replaced at the time with a temporary, lower-efficiency detector. 


s and were adjusted during the campaign as the nights grew 
shorter. Longer integrations were usually split into two indi¬ 
vidual exposures to facilitate cosmic-ray rejection. 

Each AGN was observed at a fixed slit position angle (PA) 
throughout the program in order to maintain consistency in the 
amount of spatially extended host-galaxy emission within the 
slit (TableQ}. Additionally, the flux-scaling method applied to 
the spectra relies on the assumption of a constant [O III] flux, 
and observing at a fixed PA ensures consistency in the por¬ 
tion of the narrow-line region subtended by the slit, which can 
be important for low-redshift AGNs having spatially extended 
emission. The PA for each object was chosen to be close to 
the parallactic angle (iFilippenkdl 1982 1 during the portion of 
the campaign when the object would be observed at relatively 
higher airmass. Since the observations were not carried out 
with the slit at the parallactic angle, wavelength-dependent 
slit losses have an impact on the quality of the relative flux 
calibration, although those losses are somewhat mitigated by 
the use of the 4"-wide slit. Slit losses contribute spurious scat¬ 
ter to the continuum light curves measured from the blue end 
of the spectra and to the Ha light curves measured from the 
red-side data. For Zw 229-015, the slit PA was set to 56.°5 in 
order to include a foreground star in the slit_and for consis¬ 
tency with previous measurements (iBarth et al .11201 lab . 

On clear nights we attempted to observe each AGN in the 
sample, but Zw 229-015 (at right ascension 19 hours) was 
only occasionally observed during the first month of the pro¬ 
gram. As the campaign progressed and nights became shorter, 
exposure times were shortened for most objects. We also 
ended monitoring for Mrk 141, Mrk 486, Mrk 817, and Mrk 
841 about two weeks before the end of the campaign in order 
to conserve observing time for higher-priority targets. The 
weather at Lick was worse than average during the spring of 
2011. Of the 69 nights, 20 were entirely lost to bad weather. 
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and another 12 nights suffered > 50% losses of observing 
time to clouds or high humidity. Only 23 nights were fully 
usable with no significant weather-related losses, while the 
remaining 14 nights were mildly to moderately impacted by 
clouds. Consequently, the number of observations per AGN 
was relatively small, and for most targets we obtained be¬ 
tween 30 and 40 observations during the 69-night program. 
Table □ lists the number of exposures, median airmass, and 
median S/N at 4600-4700 A for each object. During partially 
usable nights, the objects observed at highest priority included 
Mrk 40, Mrk 50, Mrk 1511, and NGC 4593. 

In order to extend our light curves for two AGNs, we also 
requested additional observations from other observers using 
the Kast spectrograph: Mrk 50 from January through March 
of 2011, and Zw 229-015 in June and July. The Mrk 50 data 
(including obser vatio ns by other observers) were previously 
presented by iBarth et al.l (1201 lbi) . Unfortunately, the 12 ob¬ 
servations taken prior to the start of our campaign suffered 
from large sampling gaps owing to winter weather, and do 
not add much additional information to the lag measurements. 
ForZw 229-015, three additional observations were taken 20- 
23 days after the end of our main campaign. 

During each night, time was also set aside for one spec¬ 
troscopic observation of a bright supernova or other tran¬ 
sient. Observations taken during our program have been in¬ 
cluded in publications exa mining the s pectral evolution o f SN 
20100 (ISmith et al.l2012 ). SN 201 lav (Foley et al.iCOTl) . SN 
201 lby (ISilverma n et al.ll 2013l). and the young stellar object 
PTF lOnvg dHillenbrand et al.S2013l) . and have been described 
in Central Bureau Electronic Telegrams (CBETs) 2681 (SN 
201 lay), 2699 (SN 201 lbp), 2701 (SN 201 lbq), 2702 (SN 
201 lbr), 2712 (SN 201 lcc), 2716 (SN 201 lcf), 2721 (SN 
201 lcj), and 2752 (SN 201 ldt)0 

4. DATA PROCESSING 
4. 1. Spectroscopic Data Reduction 

The initial data reduction followed standard procedures in¬ 
cluding overscan subtraction, two-dimensional bias subtrac¬ 
tion (for the red-side CCD only), flat-fielding, and cosmic-ray 
cleaning using the LA-Cosmic routine (Ivan Dokkumil200ll) 
implemented as an IR AI FI script. The IRAF response task 
was used to produce normalized flat-field frames having pixel 
values close to unity. Spectral extractions were done with 
the IRAF apall task. An unweighted extraction of width 
10"3 was adopted in order to minimize the impact of nightly 
seeing variations. The wavelength scale was established by 
applying a polynomial fit to the line-lamp spectra observed 
each afternoon, and then applied to the extracted spectra. For 
each wavelength-calibrated spectrum, a small final linear shift 
to the wavelength scale was applied based on offsets mea¬ 
sured from strong emission lines in the corresponding night- 
sky spectrum. Flux calibration and removal of tell uric ab - 
sorption features (following the procedure of lWade & Hornel 
19881) were carried ou t using IDL routines as described by 
Matheson et al. ( 2000 ). The calibrated spectra were rebinned 
to a linear wavelength scale of 1.0 A pixel" 1 on the blue side 
and 2.0 A pixel -1 on the red side. When two exposures of an 

29 CBETs are distributed via 

http://www.cbat.eps.harvard.edu/index.html . 

30 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As¬ 
tronomy, Inc., under cooperative agreement with the National Science Foun¬ 
dation. 


AGN were taken on the same night, they were combined to 
produce a final weighted average spectrum. 

Error spectra produced by the apall task were propa¬ 
gated through all steps of the calibration pipeline. For spectra 
of these spatially extended AGNs, it is necessary to use the 
unweighted extractions because optimal extractions ( Horne 
11986 ) have a tendency to truncate the peaks of strong emission 
lines such as [O III] when they have a different spatial extent 
than the surrounding continuum flux. However, the apall 
error spectrum actually gives the pixel uncertainties (owing to 
photon counting and readout noise) on the optimally weighted 
extraction, not on the unweighted extraction. For our high- 
S/N observations (typically with S/N « 50-100 pixel -1 ), the 
differences between the unweighted and optimally weighted 
extractions are very small, and the error spectra will only 
slightly underestimate the actual photon-counting uncertain¬ 
ties. As described below, systematic uncertainties in flux cal¬ 
ibration are often a major or dominant source of noise in the 
spectroscopic light curves. This additional source of error is 
not accounted for in the extracted error spectra and must be 
added separately to the error budget after the spectral rescal¬ 
ing procedure is applied. 

4.2. Photometric Scaling 

Conditions at Lick are rarely photometric during the spring 
season, and the flux-calibrated spectra do not share a con¬ 
sistent overall flux scale. The wavelength scales of the cal¬ 
ibrated spectra also exhibit small offsets from night to night, 
caused by random miscentering of the AGN in the spectro¬ 
graph slit. Before carrying out the emission-line measure¬ 
ments, the wavelength and flux scales of the spectra must be 
aligned. 

The blue-side data were treated as follows for each AGN in¬ 
dividually. First, wavelength shifts between the nightly spec¬ 
tra were measured using a cross-correlation procedure, sam¬ 
pling in steps of 0.1 pixel. The wavelength shifts determined 
from cross-correlation were removed from the wavelength 
vector of each spectrum, and the aligned spectra were then 
averaged to produce an initial mean reference spectrum. 

In order to place the reference spectrum on an approximate 
absolute photometric scale, we examined the nightly observ¬ 
ing logs and selected nights listed by the observers as being 
either photometric or very clear; these candidate clear nights 
numbered only 14 during the 69-night campaign. For each 
AGN, the [O III] A5007 flux was measured from each candi¬ 
date clear night’s spectrum by integrating the line flux above 
a straight-line fit to adjacent continuum regions on either side 
of the line. Then, the standard deviation of these line fluxes 
was calculated, and > 2 a outliers were discarded from the list 
of candidate photometric observations. This process was then 
repeated until there were no remaining 2o outliers in [O III] 
flux, and the mean [O III] flux from the remaining spectra was 
taken to be the best estimate of the true line flux. The num¬ 
ber of final “clear-night” spectra varied for each AGN and 
ranged between 6 and 9. Table [2] lists the [O III] fluxes based 
on this procedure; the quoted uncertainty in the table is the 
standard deviation of the fluxes measured on candidate pho¬ 
tometric nights after excluding outliers. The reference spec¬ 
trum was then multiplied by a scaling factor to normalize its 
[O III] flux to this photometric value. We caution that these 
[O III] fluxes should be considered only as estimates of the 
true fluxes, given the scarcity of truly photometric conditions 
at Lick. 

Then, the procedure described by 
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Observed Wavelength (A.) 


Fig. 1.— Mean spectra of the sample including data from the blue-side and red- 

Ivan Groningen & Wanders ( 199 2.) was applied to the data to 
align the wavelength and flux scales under the assumption that 
the [O III] flux is intrinsically constant during the duration of 
the monitoring campaign. Long-duration monitoring of NGC 
5548 confirms that this is a reasonable assumption, although 
[O III] flu x variati ons have been det ecte d on timescales of 
years dPeterson et al.1120131) . The Ivan Groningen & Wandersl 
(1992) procedure applies a scaling factor, a linear wavelength 
shift, and a Gaussian broadening to the individual night’s 
spectrum, and finds the parameter values which minimize 
the value of x 2 when fitting a low-order polynomial to the 
difference between the reference spectrum and the adjusted 
individual spectrum, over a small wavelength range including 
[O III] and some redward continuum. Since the spectral 
focus and line-spread function can vary from night to night, 
and either the individual spectrum or the reference spectrum 
might have a narrower line-spread function, the procedure 


-side cameras. The displayed wavelength range is truncated at 7500 A. 

is repeated with the Gaussian broadening kernel applied to 
the reference spectrum instead of the individual spectrum, 
and the final scaled spectrum is taken to be the better of the 
two options in terms of the x 2 value obtained. Instead of 
using a grid search to find the best-fit_values of the three 
free parameters as described by Ivan Groningen & Wandersl 
(fl992h . our method uses downhill simplex optimization, 
implemented with the IDL amoeba procedure, to minimize 
X 2 - We found that over the full time series of spectra for a 
given AGN, the shifts in wavelength scale needed to align the 
spectra were at the level of ~ 1 A rms. The final blue-side 
spectra are rebinned to a uniform wavelength scale of 1.0 A 
pixel -1 with integer wavelength values in A. 

To test the effectiveness of the spectral scaling, we can ex¬ 
amine the scatter in [O III] flux values for each AGN over the 
course of the campaign. Since the scaling procedure essen- 
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TABLE 2 

[O m] and Featureless Continuum Fluxes 


Object 

/([O III]) 

o-„ x ([0 III]) 

fx (5100 A) 

Mrk 40 

68.8±6.2 

0.007 

0.56 ±0.08 

Mrk 50 

16.3 ±2.1 

0.033 

1.20 ±0.20 

Mrk 141 

34.1 ±3.7 

0.016 

1.09 ±0.16 

Mrk 279 

157.6 ± 7.7 

0.006 

0.93 ±0.10 

Mrk 486 

65.4±4.3 

0.019 

2.15 ±0.26 

Mrk 493 

25.8 ±1.6 

0.023 

1.43 ±0.17 

Mrk 504 

17.3 ±0.5 

0.015 

0.57 ±0.06 

Mrk 704 

114.7 ± 14.5 

0.005 

2.32 ±0.37 

Mrk 817 

143.7 ±18.7 

0.013 

5.69 ±0.92 

Mrk 841 

270.4± 19.1 

0.009 

3.63 ±0.45 

Mrk 1392 

180.0 ±28.8 

0.006 

2.41 ±0.45 

Mrk 1511 

31.3 ±0.5 

0.015 

0.86 ±0.09 

NGC 4593 

144.4 ±5.8 

0.023 

2.64 ±0.29 

PG 1310-108 

128.2 ±2.7 

0.007 

1.66 ±0.17 

Zw 229-015 

37.4 ±1.3 

0.007 

0.56 ±0.06 


Note. — [O III] line fluxes are given in the observed 
frame in units of 10 -15 erg cm -2 s -1 , and uncertainties rep¬ 
resent the standard deviation of [O III] fluxes measured on 
candidate photometric nights. The quantity cr nx is the nor¬ 
malized excess scatter in the [O III] light curve, as defined in 
E3 For Mrk 50, data points taken before the start of the 
main campaign were excluded from the computation of cr nx , 
because the earlier observations were taken with heteroge¬ 
neous spectrograph setups and had a much larger flux scaling 
scatter. The AGN featureless continuum flux density is given 
in the AGN rest frame at 5100 A in units of 10~ 15 erg cm -2 
s' 1 A -1 , and is measured from the power-law component of 
the mean spectrum decomposition and corrected for Galactic 
extinction. 

tially matches the [O III] profile for each nightly spectrum to 
the reference spectrum, but does not explicitly force all [O III] 
fluxes to be identical, the dispersion among [O III] fluxes for a 
given AGN gives a measure of the final random uncertainties 
in the flux calibration after the scaling procedure has been ap¬ 
plied. We use the normalized excess variance to quantify 
the light-curve variance over and above the amount expected 
from the propagated measurement uncertainties: 

1 N 

= (i) 

where N is the number of observations in the light curve, (/) 
is the mean flux, and /j and d; are the individual flux mea¬ 
surements and their uncertainties from photon counting and 
detector read noiseFI Then, cr nx is a measure of the frac¬ 
tional rms scatter in the light curve over and above the amount 
expected from the propagated uncertainties; we refer to this 
quantity as the “normalized excess scatter.” We measured er nx 
from the [O III] light curves for each AGN, and the values 
are listed in Table [2] The values of er nx range from 0.005 to 
0.033 (i.e., the excess scatter in the light curves ranges from 
0.5% to 3.3% of the mean [O III] fluxes). The objects with the 
highest [O III] photometric scatter tend to be AGNs having 
low [O III] e quiva lent width, such as Mrk 50 and Mrk 493, for 
which the Ivan Groningen & Wanders! d 19921) scaling method 
does not work optimally. 

We interpret cr nx as a measure of the level of random error 
in the final flux normalization of the spectra. For measure- 

31 For clarity, throughout this paper we use the letter S to denote measure- 
ment uncertainties, and a to denote dispersions of a sample of measurements 
or dispersions (widths) of line profiles. 


ments of reverberation lags or other quantities from the light 
curve, this source of error can be incorporated into the mea¬ 
surement error budget by adding tr nx x (/) in quadrature to 
the photon-counting uncertainties on each light-curve point 
/. Since the spectral scaling is optimized for the wavelength 
of the [O III] A5007 line, the relative scaling at much longer or 
shorter wavelengths is expected to be worse than the amount 
indicated by cr nx , as a result of atmospheric dispersion and 
wavelength-dependent slit losses. 

After applying the scaling procedure, mean and rms spectra 
are constructed following iPeterson et al.l ( 2004 ). The mean 
spectrum is a simple average of the individual nightly spectra. 
The rms spectrum is calculated by taking the standard devia¬ 
tion of flux values at each pixel over the time series, and gives 
an illustration of the relative variability amplitude of different 
spectral regions. Other pres cripti ons for calculating rms spec¬ 
tra have been explored bv iPark et al.l (12012 a:). Since the rms 
spectrum includes the contribution of photon-counting noise 
and other sources of noise or random error in the data, in ad¬ 
dition to genuine AGN variability, the amplitude of the rms 
spectrum is higher than it would be for the ideal case of noise- 
free data. Simulations presented in Appendix ICl demonstrate 
that photon-counting noise will bias the widths of broad lines 
in rms spectra to values that are typically several percent lower 
than they would be in the absence of noise. 

The red-side data were generally flux-calibrated using 
standard-star observations obtained simultaneously with the 
blue-side standard-star observations. In the reduced data, the 
region of wavelength overlap between the blue-side and red- 
side spectra is 5420-5500 A. For each AGN observation, we 
normalized the red-side flux scale by applying a scaling factor 
to match it to the scaled blue-side spectrum over this wave¬ 
length range, and then produced mean and rms spectra from 
these rescaled red-side data. FigureHJdisplays the mean spec¬ 
tra of each AGN from the blue and red camera data combined. 

For a few of the AGNs, we found that Ha light curves mea¬ 
sured from these rescaled spectra had noticeably larger scatter 
in flux values than the Hu light curves. This large excess scat¬ 
ter for the red-side data probably arises from a combination of 
the nonparallactic slit orientation, the seeing-dependent ratio 
of AGN to host-galaxy continuum flux, and the fact that the 
spectral normalization is optimized for wavelengths close to 
[O III]. Large excess scatter was particularly noticeable in the 
Ha light curves of Mrk 40 and_Mrk 279. For both of these ob¬ 
jects, we applied the Ivan Groningen & Wanders! 11992!) scal¬ 
ing method separately to the red-side spectra, using the [O I] 
A6300 line as a reference. This significantly improved the 
quality of the Ha light curves for these two objects, although 
a few extreme outlier points were still present in the data. For 
other objects, we found that either rescaling using [O I] did 
not noticeably improve the Ha light curves, or, in some cases 
such as Mrk 50, [OI] was too weak to use as a reference line 
in the spectral scaling. We experimented with using the [S II] 
doublet as the reference line, but this produced poor results 
because [S II] is very weak in most of these AGNs and is usu¬ 
ally blended with the red wing of Ha. 

4.3. Spectral Fitting 

Traditionally, measurements of emission-line fluxes for re¬ 
verberation mapping have been carried out by selecting two 
regions nearly free of emission lines on either side of the 
emission line of interest, then fitting and subtracting a straight 
line to model and remove the continuum underlying the emis- 
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sion line. While this procedure is simple to implement and 
can generally produce accurate light curves for strong emis¬ 
sion lines, there are several drawbacks associated with it. A 
straight-line fit does not adequately represent the shape of the 
starlight continuum, which often makes a significant contribu¬ 
tion to the spectra of Seyfert 1 nuclei (particularly when ob¬ 
served through a wide slit). Stellar absorption features in the 
host-galaxy starlight produce wiggles and bumps that become 
superposed on the broad H/3 line profile (including stellar H/3 
absorption), and these contaminating features are not removed 
when a straight-line continuum is assumed. Furthermore, 
while there are some regions around H/3 where the flux from 
other emission lines is relatively low, there may still be con¬ 
tributions from weak features such as Fe II or He II in the se¬ 
lected continuum fitting regions, and it is effectively impossi¬ 
ble to define a true continuum region that is entirely devoid of 
AGN emission features. These issues can compromise mea¬ 
surements of emission-line light curves and broad-line widths 
or profile shapes in the mean and rms spectra. Measurements 
of velocity-resolved lags in the H/3 line can be compromised 
by blending with He II, He I, and Fe II if these features are 
not removed before measurement of the reverberation signals, 
and biased determination of the true continuum level owing to 
contamination by blended lines or starlight could affect mea¬ 
surements of the lag in the faint high-velocity wings of H/3. 
In addition to strong lines such as H/3, it is also of interest 
to measure light curves for weaker features including Fe II or 
He II, but these lines can be severely blended with each other 
as well as with starlight and other spectral components, and a 
straight-line continuum fit is usually not adequate for isolating 
the contributions of weaker broad lines to the spectrum. 

In order to address these problems, we have implemented a 
spectral fitting method which, when applied to each nightly 
spectrum from a reverberation campaign, decomposes the 
data into independent spectral components. Fitting of AGN 
spectra to deblend line and continuum components (including 
the broad Fe II blends) has a long history with many differ- 


ent implementations in the literature (e. 

g., IWills et al.l 

1985; 

Dietrich et al.1 2002); |Woo et al. 2006; 

Greene & Ho 

2007; 

Shen et all 2008; Kovacevic et all 201C 

), but only recently 


have such methods been applied to time-series spectroscopic 
data from reverberation campaigns. The procedure described 
below is an updated version of the met hod applied to our data 
in earlier papers in this series dBarth et al.li201 lbl 120131) . with 
several minor modifications and improvements. The method 
described bv lPark et all <*2 012a ) uses a different fitting proce¬ 
dure and a different set of model components, but we find that 
measurements of H/3 light curves and other parameters are 
very similar when carried out using the results of the two in¬ 
dependent fitting codes when applied to the same set of input 
spectra. 

The primary goal of our procedure is to provide an accurate 
fit to the region surrounding H/3 so that the broad H/3 profile 
can be isolated from other features. Consequently, accurate 
fitting over a small wavelength region is a higher priority than 
fitting the entire available wavelength range of th e observa- 
tions. Fits to a larger wavelength region (e.g., [Denney et all 
2009 a) would include several additional emission lines that 
are of little direct interest for the reverberation measurement, 
and extending the fitting region blueward could compromise 
the accuracy of the fit in the H/3 region owing to the com¬ 
plexity of the spectrum at bluer wavelengths. We fit only 
the blue-side spectra, and the redward limit of the fits is set 
by the red end of the blue-side CCD data. Fits were car¬ 


ried out over a rest wavelength range of approximately 4200- 
5300 A, with slight differences for each AGN depending on 
redshift, and the fitting procedure is app lied to the data after 
scaling by the Ivan Groningen & Wandersl (119921) procedure. 
The Levenberg-Marquardt routines in the IDL mpfit pack¬ 
age (lMarkwardtll2009l) are employed for \ 2 minimization to 
optimize the fits, using the propagated error spectra for deter¬ 
mination of x 2 - The model components and parameters are as 
follows. 

Starlight: W e used an 11 Gyr, so l ar met allicity, single-burst 
spectrum from iBruzual & Charlol ( 2003 ) to model the host- 
galaxy starlight. We additionally carried out tests in which 
a younger (290 Myr) population component was added. Its 
contribution was generally small and not well constrained, 
and it was omitted in the final fits. During the fitting pro¬ 
cess, the starlight spectrum was rebinned to a log-linear wave- 
lenth scale, broadened by convolution with a Gaussian ker¬ 
nel, and then rebinned again to match the linear wavelength 
scale of the data. Free parameters include the flux normal¬ 
ization of the starlight spectrum, the redshift, and the veloc¬ 
ity width of the broadening kernel. S ince the spectral resolu¬ 
tion of the IBruzual & Chariot] (1 20031) spectra differs from the 
Lick data, the convolution kernel width does not directly cor¬ 
respond to the actual stellar velocity dispersion in the host 
galaxy. The velocity-broadening dispersion parameter was 
limited to a maximum value of 350 km s -1 in order to prevent 
it from increasing to unphysically large values; this limit was 
reached only for Mrk 486, in which stellar absorption features 
are essentially invisible over the wavelength range of the fits. 

AGN featureless continuum: The AGN continuum was 
modeled as a power law, with free parameters including the 
flux normalization at rest wavelength 5100 A and the power- 
law spectral index. Over the wavelength range that was fit¬ 
ted, Balmer continuum emission does not contribute to the 
spectrum. We do not ascribe any physical meaning to night- 
to-night changes in the best-fitting spectral index, since these 
can be caused by wavelength-dependent slit losses resulting 
from differential atmospheric refraction and variable seeing. 

[O 1111: A fourth-order Gauss-Hermite function 
(Ivan der Marel & Franxlll993 ) was used to model the [O III] 
A5007 profile. Free parameters included the amplitude, 
centroid wavelength, line dispersion, and the higher-order 
moments hj and /z 4 . The 4959 A line was modeled with the 
same velocity profile and an amplitude fixed to 1/3 of that the 
5007 A line, and a fixed rest-frame wavelength separation of 
47.9 A. A fourth-order Gauss-Hermite function proved to be 
adequate for modeling the [O III] profiles well in all but one 
of the AGNs; as described below, we modified the set of fit 
components for Mrk 279. 

H/3: For the broad component of H/3, a fourth-order Gauss- 
Hermite function was used. This provides sufficient flexibility 
to fit the broad-line profiles for nearly all of the AGNs well, 
with the exception of objects such as Mrk 141 having distinct 
shoulders or humps on their profiles. Free parameters include 
the broad component centroid, dispersion, I 13 , /z 4 , and ampli¬ 
tude. The narrow component was set to have a profile identi¬ 
cal to [O III], but a small shift was allowed in the wavelength 
centroid of the narrow component as another free parameter. 

He II: The He II A4686 line was modeled as the sum of 
two Gaussians representing the broad and narrow compo¬ 
nents. For each component, free parameters included the cen¬ 
troid, dispersion, and amplitude. Initial trial fits indicated 
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that the broad and narrow components sometimes had offset 
centroid wavelengths, and for some AGNs the fits failed no¬ 
ticeably when the two components were forced to have the 
same centroid. In some cases this might be the result of broad 
He I A4713 being present in the data but not included as a 
fit component. If the broad He II model component incor¬ 
porated some He I A4713 flux then this would induce a spu- 
r ious but sm all redward shift relative to its true wavelength. 
IVeron et~akl d 2002 l) illustrate examples of fits that include both 
He II A4686 and He I A4713, but since these features are badly 
degenerate in broad-lined AGNs we do not attempt to include 
the He I A4713 line in our fits as a separate component. 

He I: We included the possible contributions of He I broad 
lines at 4471, 4922, and 5016 A rest wavelength. In order to 
keep the number of free parameters to a minimum for these 
relatively weak lines, we used Gaussian profiles and followed 
IVestergaard & Petersonl (120051) in setting the 4922 and 5016 A 
lines to have the same amplitude, but we did not force any re¬ 
lationship between the amplitudes of these lines and the 4471 
A line. All three He I lines were assumed to have identical 
velocity widths, and their velocity centroids were set to be 
identical to that of H/3. The He I components thus required 
only three additional free parameters to represent the ampli¬ 
tudes and widths of these features. The He I 4922 and 5016 
A lines are in the “red shelf” of H/i, and are blended with 
the red wing of H/3, the [O III] lines, and with Fe II features 
(IVeron et al .1 120021) . In fact the 4922 and 5016 A lines are 
nearly degenerate with two Fe II features that appear in the 
iron templates. The amplitude of these two features in the fits 
depends strongly on which Fe II template is chosen. While the 
sum of the Fe II and He I fluxes in this spectral region is fairly 
well determined by the fits, we do not attempt to extract He I 
light curves because of the ambiguity in separating these fea¬ 
tures from the Fe II emission. For most AGNs, the best-fitting 
amplitude of the He I A4471 line is either zero or negligibly 
small, regardless of Fe II template choice. 

Fe II: The blends of broad Fe II emission were mod- 
eled using templa t e spec tr a, including those describ ed by 
Boroson & Greenl (119921) . IVeron-Cettv et alJ (120041) . and 
Kovacevic et al.l (120101) . We carried out fits separately with 
each of these templates in order to examine systematic differ¬ 
ences in the results. Template spectra were convolved with 
a Gaussian kernel in velocity space, similar to the broaden¬ 
ing of the starlight spectrum described above. In each case, 
free parameters included the Fe II velocity shift relative to 
broad H/3, the broadening kernel width, and the normali zatio n 
of the broa dened template spe c trum. The iBoroson & Greenl 
(11992b and IVeron-Cettv et al.l (120041) templates are mono- 
lithic and only one nor malization paramet er is required. The 
Kovacevic et al.l (120101) template (see also iShanovalova et alJ 
20121 ) is composed of five separate templates representing dif¬ 
ferent emission multiplet groups, and a normalization param¬ 
eter i s requ ired for each one. In a given spectral fit using the 
iKovacevic et al.l (120101) templates, each of the five templates 
was assumed to have the same velocity broadening and veloc¬ 
ity shift. The lKovacevic et aD (1201 01) templates are distributed 
with a discrete sequence of velocity broadenings, and we used 
the narrowest templates (corresponding to 700 km s -1 broad¬ 
ening) and then applied a Gaussian velocity broadening dur¬ 
ing the fitting process as described above. Examination of the 
fits for a very narrow-lined AGN such as Mrk 493 clearly il¬ 
lustrate the contribution of Fe II in the red wings of both H/3 
and [O III] A5007, while in more typical broad-lined AGNs 



Fig. 2.— Illust ration of the differences between the F e II templates of 
Boroson & Greenl < 19921) (blue), IVeron-Cettv et al] <20041) (magenta), and 
Kovacevic et alJ <20101) (red) when fitted to the spectrum of Mrk 493. The 
Mrk 493 mean spectrum, after subtraction of the AGN featureless continuum 
and starlight components, is shown in black, and has been shifted upward by 
0.7 units for clarity. The colored curves show the Fe II templates with veloc¬ 
ity broadening and normalization as determined from the best fit to the mean 
spectrum. 

these Fe II features are more difficult to identify individually 
and their contribution to this spectral region may be less ob¬ 
vious unless a multicomponent fit is carried out. 

Reddening : All model components were reddened by a 
ICardelli et all (119891) reddening law, with E(B-V) as a free 
parameter in the fit. This parameter accounts for possible red¬ 
dening within the AGN host galaxy as well as Galactic red¬ 
dening, and also accounts for possible errors in relative flux 
calibration resulting from the nonparallactic slit orientation or 
other wavelength-dependent flux calibration errors. For some 
AGNs, the best-fitting values of E(B-V) varied significantly 
from night to night, and we do not interpret the fitted parame¬ 
ter values as actual measurements of the reddening toward the 
AGNs. In principle, the AGN continuum and BLR, NLR, and 
host-galaxy light could have different line-of-sight reddening, 
but the data would not provide sufficient constraints to obtain 
meaningful values of reddening for these components inde¬ 
pendently. 

In the H/3 and [O III] models, the wings of Gauss-Hermite 
functions can extend to negative values for certain values of 
lij and h. 4 . Negative values were clipped and replaced with 
zero. This only affects the distant wings of line profiles and 
has a very minor effect on the fits. 

Other emission lines are present in the spectrum that could 
be fitted by adding components and free parameters to the 
model, but we chose to restrict the number of components 
to those listed above in order to avoid either slowing down 
the fitting procedure or introducing additional fitting degen¬ 
eracies. Instead, other emission features were accounted for 
simply by setting a wavelength region around those features 
to have zero weight in the fits. The largest feature that was 
masked out was the blend of H 7 and [O III] A4363. This 
blend typically has a complex shape that, if included in the 
fit, would require three overlapping Gauss-Hermite or Gaus¬ 
sian components and potentially nine or more additional free 
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Fig. 3.— Fits to the blue-side mean spectra. The observed mean spectra are shown i n black. Model components include starlight (purple), the AGN featureless 
continuum (green), H/3 (magenta). He II (cyan). He I (orange). IKovacevic et all <20101) Fe II template (grey), and [O III] (blue). The sum of all model components 
is shown in red, superposed on the observed spectrum. The residual spectrum (data — total model) is shown in black and is offset downward to negative values 
for clarity. This residual contains emission features including H 7 , [O III] A4363, and [NI] A5200 that are not incorporated as components in the model. 


parameters. The masked region corresponded to 4280-4400 
A (rest wavelength) but was adjusted in a few cases for objects 
having very broad lines. We also masked out a small region 
around the narrow [NI] lines at 5199 and 5201 A. 

Our procedure fits all of these spectral components to the 
data simultaneously, in contrast to other methods which fit 
and remove continuum and Fe II components prior to fit - 
ting the H/3 profile (e.g JShen et al. 2008 ; [Park etaTlf2(rf2a . 
Fits were carried out separately using each of the three Fe II 


templates. The total number of free parameters in the fits 
is 29 when the monolithic Fe II templates are used, or 33 
when the multicomponent IKovacevic et all (i2010l) templates 
are used. For each AGN, the mean spectrum was fitted first, 
and then the best-fit parameters determined for the high S/N 
mean spectrum were used as the initial parameters for the 
fits to each nightly spectrum. In most cases, fits using the 
IKovacevic et al.l i2010l) templates returned the best xl values, 
and we use these templates for the final fit results presented 
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Fig. 4.— Examples of fits to the H/3 line in single-epoch observations. The 
displayed data are from observations taken on 2011-04-27. In each panel, 
the decomposed H/3 profile obtained by subtracting all model components 
other than H/3 from the spectrum is shown in black. This includes both the 
narrow and broad components of H/3. The broad component spectrum after 
additionally subtracting the narrow H/3 model component is shown in grey. 
Blue and red curves show the best-fitting broad H/3 and total H/3 (broad + 
narrow) model components. 

in this paper. Figure [2] illustrates the systematic differences 
between the three templates when fitted to the spectrum of 
Mrk 493, which has the strongest and narrowest Fe II lines of 
the AGNs in our sample. The template differences are most 
pronounced at wavelengths underlying the broad H 7 and H/3 
lines and the H/3 red shelf region. 

The fitting procedure was modified for Mrk 279. Its [O III] 
lines are unusually broad with strong asymmetric blue wings 
which contribute substantially to the “red shelf” region red- 
ward of H/3, and the 4th-order Gauss-Hermite model gave a 
poor fit to the [O III] doublet leaving very strong residuals in 
this region. For this object, we removed the He I components 
from the model (since they did not appear to be necessary 
in this case) and replaced them with two additional Gaus¬ 
sian components representing the blue wings of the [O III] 
AA4959,5007 lines. The H/3 model also required modifica¬ 
tions. In order to prevent the broad H/3 width from blowing 
up to arbitrarily high values due to fitting degeneracy in the 
red wing, we restricted the broad component model to a < 30 
A. We also found that the narrow component of H/3 had a 
significantly different profile from [O III] in Mrk 279, and al¬ 
lowed the width of narrow H/3 to vary rather than tying it to 
the [O III] width. This resulted in cr(H/3 n )« 0.75cr([O III]). 

Th e best fits to the mean spectrum of each AGN, using the 
iKovacevic et alJ (120101) template, are shown in Figure [3] 

For each spectrum, the best-fitting model components are 
saved. Additionally, we save the decomposed H/3 profile, 
made by subtracting all model components other than H/3 
from the original spectrum. This spectrum also contains the 
H 7+[0 III] blend as part of the fitting residual. We save two 
versions of this spectrum, one including and one excluding 
the narrow H/3 component. These decomposed spectra (rather 


than the noise-free models) are used for measurement of the 
H/3 light curves and line widths. Figure[4]illustrates examples 
of the decomposed H/3 profiles and component models for six 
AGNs from a single night’s data. 

The fitting procedure returns estimates of uncertainties on 
each of the free parameters. However, for measurements of 
emission-line light curves, what is needed is an estimate of 
the uncertainty on the flux density at each pixel for a partic¬ 
ular spectral component such as H/3. We employed a Monte 
Carlo procedure to estimate the flux uncertainties. For each 
spectrum, after finding the best model fit, we created 100 
modified versions of the spectrum by adjusting each pixel by 
a random Gaussian deviate drawn from a distribution with a 
Gaussian dispersion equal to the amplitude of the propagated 
error spectrum at that pixel. The spectral fit was then repeated 
for each of the 100 randomly deviated spectra, and the in¬ 
dividual fit components were saved. For each spectral com¬ 
ponent, the value of the error spectrum at a given pixel was 
taken to be the standard deviation of the flux values from the 
100 re-fitted versions of that spectral component. Similarly, 
we produced error spectra for the decomposed H/3 profiles by 
taking the standard deviation of the 100 Monte Carlo realiza¬ 
tions of the decomposed profile (i.e., the deviated spectrum 
minus all best-fit components other than H/3). 

From this procedure, we found that the error spectrum of 
the decomposed H/3 profile has an amplitude at each spectral 
pixel that is nearly equal to that of the original error spec¬ 
trum, but noisier as a result of the limited number of Monte 
Carlo iterations. Based on the results of this procedure, we 
adopt the original error spectrum as our best estimate of the 
pixel uncertainties on the individual emission-line component 
spectra that we use to measure light curves. The error spec¬ 
tra do not incorporate possible systematic uncertainties due to 
differences between Fe II templates or the choice of a partic¬ 
ular host-galaxy starlight model. 

Using the fit results, we constructed revised mean and rms 
spectra. Figure 0 shows the original mean and rms spectra 
along with mean and rms spectra constructed after subtrac¬ 
tion of the best-fitting starlight and AGN featureless contin¬ 
uum components from each nightly spectrum. In some cases 
such as Mrk 279, the H/3 profile in the rms spectrum changes 
significantly after subtraction of the continuum components; 
similar results have previously been discussed by [Park et all 
(l2012af) . Additionally, for some AGNs the He II profiles in 
rms spectra are significantly altered by removal of the contin¬ 
uum components; these changes are particularly prominent in 
Mrk 504 and Mrk 1511. Another version of the rms spectrum 
was constructed by subtracting all model components other 
than broad H/3 from the nightly spectra. This version isolates 
the variability of the broad H/3 line itself, and was used to 
measure the width of the H/3 rms profile as described in Sec¬ 
tion [ 6 ] Figure [ 6 ] shows multiple versions of the rms spectrum 
for Mrk 40 to illustrate the changes when different spectral 
components are subtracted prior to calculating the rms flux. 
One advantage of subtracting narrow emission-line compo¬ 
nents prior to constructing the rms spectrum is that the resid¬ 
ual noise from these intrinsically nonvariable features is sig¬ 
nificantly reduced. This is clearly apparent for the [O III] lines 
in Mrk 40 (Figure [ 6 |. 

In future work, a variety of possible extensions and im¬ 
provements to the fitting method can be explored. These in¬ 
clude (a) addition of model components for other emission 
lines such as H 7 and [O III] A4363 instead of simply masking 
them out from the fit; (b) broadening the wavelength range 
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Fig. 5.— Mean and rms spectra for the blue-side fitting region. In each panel, the upper subpanel shows shows the standard mean spectrum (blue), and a “line- 
only” mean spectrum constructed after subtraction of the AGN featureless continuum and starlight components from each nightly spectrum (red). The lower 
subpanel shows the standard rms spectrum (blue), and an rms spectrum constructed after subtraction of the AGN featureless continuum and starlight components 
from each nightly spectrum (red). 


of the fit blueward down to the near-UV atmospheric cutoff, 
and redward to Ha; (c) inclusion of multiple-age stellar pop¬ 
ulations; (d) addition of a Balmer continuum component for 
fits extending below the Balmer break; (e) use of higher-order 
models to achieve more accurate fits to strong lines such as 
[O III] A5007 and for irregular or double-peaked broad Balmer 
lines; (f) exploration of different methods for constraining the 
decomposition in the H/3 red shelf region; (g) allowing for 
the possibility that different reddening values or even differ¬ 


ent reddening laws might apply to spectral components orig¬ 
inating at different spatial scales. The possible benefits of 
these improvements must however be weighed against the in¬ 
creased complexity that comes with adding more free param¬ 
eters and the degeneracies among them, particularly for the 
continuum components. Furthermore, since the primary goal 
is often to achieve the most accurate fit to the wavelength re¬ 
gion surrounding H/3, care must be taken to ensure that fitting 
a broader wavelength range does not worsen the fit in this 
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FlO. 6.— Different versions of the rms spectrum for Mrk 40. The top panel 
shows the standard rms spectrum constructed using the nightly total flux spec¬ 
tra. The second panel shows the rms after removing the AGN continuum and 
host-galaxy starlight components from each nightly spectrum. In the third 
panel, narrow emission lines have also been removed ([O III] AA4959.5007, 
H/3, and He II), and in the bottom panel, the broad He II and Fe II components 
are also removed from the nightly spectra. This last version of the rms spec¬ 
trum, with all model components removed other than broad H/3, is the one 
used to measure the width of broad H/3 in the rms spectrum. 

region. The available set of Fe II templates remains a funda¬ 
mental limitation and source of systematic uncertainty in the 
fits, and further work to develop more flexible sets of Fe II 
templates, both theoretical and emipirical, covering the full 
UV/optical wavelength range, would be extremely valuable 
for future reverberation mapping. 

5. EMISSION-LINE AND CONTINUUM LIGHT CURVES 

Light curves for broad emission lines were measured by di¬ 
rect integration of spectra. The epoch of each observation is 
given by the Heliocentric Julian Date (HJD) determined from 
the UT date saved in the FITS image headers, using the IRAF 
set jd routine. To measure the H/3 and H 7 light curves for 
each night’s observation of an AGN, we began with the ob¬ 
served spectrum and subtracted all of the best-fitting model 
components except for the broad and narrow H/3 models. This 
left a residual spectrum with a continuum level of zero (mod¬ 
ulated by the fitting residuals) and the H/3 line and H 7+[0 III] 
blends as the remaining emission features. The flux in these 
features was integrated over wavelength ranges that were se¬ 
lected for each AGN to encompass the full width of the line 
profiles. Wavelength ranges for each light curve are listed in 
Table[3] We measure the total (broad + narrow) flux for H/3 in 
order to avoid introducing possible errors into the light curves 
due to degeneracy in decomposing the two components, con¬ 
sistent with met hods typically used for reverberation mea¬ 
surements (e.g., iKaspLetal] EoQS iBentzeTaDHoolbi). The 
narrow component adds a constant flux pedestal to the light 
curves. The H 7 measurement includes the flux of the broad 
and narrow H 7 lines as well as [O III] A4363. 

We also measured light curves for He II A4686 and for the 
Fe II blends integrated over the rest wavelength range 4400- 


4700 A. The light curves for these weak features were mea¬ 
sured from the best-fitting model component spectra for each 
night’s data, and for He II the light curves include both the 
broad and narrow components. The Fe II emission was in¬ 
tegrated over the rest wavelength ranges of 4400-4700 and 
5100-5400 A, avoiding regions in which Fe II emission is de¬ 
generate with H/3 or the possible He I emission in the H/3 red 
shelf. 

For the Ha and HA lines, which lie outside the wavelength 
range of the spectral fits, we measured fluxes from the scaled 
spectra using the traditional approach of fitting a local linear 
continuum to surrounding line-free regions on either side of 
the line. The Ha light curves include the fluxes of the narrow 
[Nil] A6548,6583lines. 

Light curves are displayed in Figures [7UI41 and the light 
curve data are listed in Tabled 

Light curves were also measured for the AGN continuum 
in the U -band spectral region, since the AGN featureless con¬ 
tinuum tends to dominate over starlight in this region and is 
more highly variable than at redder wavelengths. The mean 
continuum flux density was measured over A rest = 3500-3600 
A, and we refer to this spectroscopic continuum flux as the U s 
band. For a few objects, this wavelength range was adjusted 
in order to avoid emission lines. This measurement includes 
both AGN and starlight continuum. Table [3] lists the wave¬ 
length range that was integrated for each light curve. Since 
this wavelength range is far from the wavelength of [O III], 
which was used as the flux calibration reference for spectral 
scaling, the relative accuracy of the photometric calibration is 
sometimes poor. From the appearance of the U s light curves, 
the flux uncertainties in portions of the light curves are ~ 10 % 
or occasionally worse. The scatter in the U s light curves for 
some objects (such as Mrk 50) is highest at times close to the 
start and end of the campaign, when targets were observed at 
higher airmasses, but this trend does not appear to hold for the 
entire sample. 

From each light curve, we calculated the statistics F var and 
R max to characterize the variability amplitude, and for compar¬ 
ison with previous reverberation campaigns (e.g jBentz et al.l 
l2009bl) . The quantity F v ar is similar to the normalized excess 
scatter and is defined as 


<7/-<52ns 


(/) 


( 2 ) 


where crj is the variance of the flux values in the light curve, 
Arms is the rms m easu rement uncertainty on the fluxe_s, and (/) 
is the mean flux (Rodrfguez-Pascual et alJl997tlEdelson et al.l 
12001 . R max gives the maximum variability amplitude and is 
simply the ratio of maximum to minimum fluxes in the light 
curve. To account for the contribution of residual scatter from 
the spectral scaling to the error budget, we modified the flux 
uncertainties by taking the product of excess [O III] flux scat¬ 
ter with mean flux ( a„ x x (/)) and adding that quantity in 
quadrature to each individual flux uncertainty value prior to 
calculating F var , R max , and their uncertainties. As an overall 
measure of light-curve variability, F var is preferable to R max 
because R max is particularly sensitive to outlier values at the 
highest and lowest fluxes. 

Emission-line light curves having F var >0.1 correspond to 
strong variability from which lags can generally be measured 
well. In our sample, the objects exceeding this threshold for 
H/3 variations are Mrk 40, Mrk 50, Mrk 1511, NGC 4593, 
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TABLE 3 

Light Curve Wavelength Ranges and Variability Statistics 


Galaxy Light Curve Rest Wavelengths (A) F yar R max 


Mrk 40 U s 


Mrk 50 

Ho 

H/3 

H 7 

H5 

Hell 

Fell 

U s 

Mrk 141 

Ha 

H/3 

H 7 

H<5 

Hell 

FeU 

U s 

Mrk 279 

H/3 

Us 

Mrk 486 

Ha 

H/3 

H 7 

H<5 

Hell 

Us 

Mrk 493 

H/3 

Us 

Mrk 504 

H/3 

Us 

Mrk 704 

H/3 

Us 

Mrk 817 

Ha 

H/3 

H 7 

H<5 

Hell 

Us 

Mrk 841 

H/3 

Us 

Mrk 1392 

H/3 

Us 

Mrk 1511 

H/3 

Us 

NGC 4593 

Ha 

H/3 

H 7 

H5 

Hell 

Fell 

Us 

PG 1310-108 

Ha 

H/3 

H 7 

H<5 

Hell 

Fell 

Us 

Zw 229-015 

H/3 

Us 


Ha 

H/3 

H 7 


3500-3600 

6493-6660 

4799^1916 

4309^1387 

4074—4142 

4583^1789 

4400-4700,5100-5400 

3500-3600 

6449-6693 

4758^1944 

4270^1417 

4065^1133 

4495^1817 

4400^1700,5100-5400 
3500-3600 
4782^1925 
3500-3600 
6405-6667 
4755-4920 
4251-4416 
4047-4115 
4590^1765 
3500-3600 
4784^1928 
3500-3600 
4809^1897 
3500-3600 
4798^1914 
3500-3570 
6422-6704 
4770^1936 
4256-4411 
4052^1149 
4489^1848 
3500-3600 
4750^1944 
3500-3570 
4757-4940 
3500-3600 
4758^1932 
3500-3600 
6432-6693 
4768^1933 
4236^1430 
4062-4140 
4556^1797 

4400-4700,5100-5400 

3500-3570 

6442-6690 

4757-4926 

4242-4410 

4058^1133 

4539^1807 

4400^1700,5100-5400 
3500-3600 
4787-4932 
3500-3570 
6470-6664 
4767^1912 
4281-4388 


0.30 

±0.03 

2.78 

±0.04 

0.14 

± 0.02 

2.02 

± 0.02 

0.20 

± 0.02 

1.98 

±0.03 

0.21 

± 0.02 

2.04 

±0.04 

0.27 

±0.03 

2.61 

±0.06 

0.38 

±0.04 

3.57 

± 0.21 

0.15 

± 0.02 

2.05 

± 0.10 

0.34 

±0.03 

3.88 

±0.34 

0.10 

± 0.01 

1.53 

±0.08 

0.20 

± 0.02 

2.17 

±0.13 

0.23 

± 0.02 

2.50 

±0.17 

0.29 

±0.03 

2.93 

±0.23 

0.55 

±0.05 

12.02 ±2.09 

0.35 

±0.03 

7.28 

±1.31 

0.09 

± 0.01 

1.45 

±0.03 

0.08 

± 0.01 

1.39 

±0.03 

0.13 

± 0.02 

1.63 

± 0.02 

0.04 

± 0.01 

1.22 

± 0.01 

0.07 

± 0.01 

1.31 

± 0.01 

0.07 

± 0.01 

1.34 

± 0.02 

0.09 

± 0.01 

1.45 

±0.04 

0.30 

±0.04 

2.98 

± 0.21 

0.07 

± 0.01 

1.40 

±0.04 

0.02 

± 0.01 

1.10 

±0.03 

0.04 

± 0.01 

1.20 

±0.04 

0.03 

± 0.01 

1.17 

±0.04 

0.07 

± 0.01 

1.36 

±0.03 

0.04 

± 0.01 

1.24 

±0.03 

0.13 

± 0.02 

2.02 

±0.03 

0.07 

± 0.01 

1.37 

± 0.01 

0.05 

± 0.01 

1.23 

± 0.01 

0.08 

± 0.01 

1.46 

± 0.02 

0.10 

± 0.01 

1.63 

±0.08 

0.28 

±0.04 

3.62 

±0.13 

0.08 

± 0.01 

1.40 

±0.03 

0.04 

± 0.01 

1.16 

± 0.02 

0.14 

± 0.02 

1.92 

±0.03 

0.05 

± 0.01 

1.27 

± 0.02 

0.14 

± 0.02 

1.88 

±0.03 

0.03 

± 0.01 

1.16 

± 0.01 

0.19 

± 0.02 

2.42 

±0.06 

0.09 

± 0.01 

1.45 

±0.03 

0.12 

± 0.01 

1.43 

±0.03 

0.13 

± 0.01 

1.54 ±0.04 

0.20 

± 0.02 

2.00 

±0.07 

0.54 ±0.06 

10.80 ±1.59 

0.10 

± 0.01 

1.43 

±0.03 

0.38 

±0.04 

4.54 ±0.25 

0.16 

± 0.02 

1.81 

±0.07 

0.23 

± 0.02 

2.11 

±0.08 

0.25 

±0.03 

2.15 

±0.09 

0.43 

±0.05 

4.72 

±0.50 

0.72 

±0.08 

19.28 ±4.41 

0.16 

± 0.02 

1.71 

±0.06 

0.13 

± 0.02 

1.85 

±0.03 

0.05 

± 0.01 

1.29 

± 0.02 

0.25 

±0.03 

2.65 

±0.06 

0.13 

± 0.02 

1.54 ±0.02 

0.25 

±0.03 

2.81 

±0.07 

0.27 

±0.04 

2.82 

±0.13 


NOTE. —These measurements include the constant narrow-line component contribu¬ 
tions to the Balmer-line and He II light curves, and the constant host-galaxy contribution 
to the U s light curve. The wavelength column lists the range over which flux density was 
integrated for the light curve measurement. As described in the text, the ! 1,3, H 7 , He II, 
and Fe II light curves were measured from individual components of the decomposed 
blue-side data, while U s , Ha, and Hi5 were measured from the total-flux spectra. 






15 


TABLE 4 

Light Curve Data 


Galaxy 

Light Curve 

HJD 

/ 

Sf 

Mrk 40 

Us 

5650.75 

1.081 

0.004 



5651.75 

1.112 

0.006 



5654.80 

1.191 

0.013 



5655.73 

1.146 

0.010 



5656.76 

1.127 

0.004 



5657.77 

1.111 

0.007 



5663.74 

1.386 

0.009 



5664.78 

1.481 

0.006 



5666.76 

1.515 

0.005 



5667.75 

1.567 

0.005 



5673.73 

2.101 

0.006 



5674.82 

1.848 

0.006 


Note. — Dates are listed as HJD - 2450000. Units 
are 10 1 5 erg cm ~ 2 s 1 A -1 for U s , and 10~ 15 erg cm ~ 2 
s 1 for emission-line light curves. The listed uncer¬ 
tainties 8f correspond to the propagated uncertainties 
from photon counting and detector readout noise, and 
do not include the additional contribution from resid¬ 
ual flux-scaling errors. This Table is published in its 
entirety in the electronic edition of The Astrophysical 
Journal Supplement Series. A portion is shown here 
for guidance regarding its form and content. 

and Zw 229-015. For these AGNs, we plot light curves of the 
U s flux density, the Balmer lines, and also He II and Fe II in 
those cases where there is some distinct variability above the 
level of the noise for these components. We also plot multi- 
line light curves for Mrk 279 and Mrk 704. Mrk 279 has 
a lower H/3 variability amplitude of F var = 0.07, but its light 
curve exhibits distinct up-and-down variations which provide 
useful diagnostics of the emission-line lags. Mrk 704 exhibits 
a downward trend in the continuum and emission-line light 
curves, but the nearly monotonic trend contains little short- 
timescale structure and is not optimal for determination of 
lags. For the remaining objects in the sample, all of which 
have F var (H/3)< 0.1, light curves are shown only for U s and 
H/3 in Figure [14] Aside from H/3, the light curves for other 
emission lines in these objects are too noisy to be of much 
use. A few of these AGNs, such as Mrk 141, Mrk 493, and PG 
1310-108, show distinct features in their light curves which 
may provide some limited information on the H/3 lag. Table 
[3 lists the light-curve statistics F var and R mdx for each light 
curve displayed in Figures [7HT41 For completeness, we list 
statistics measured for the U s continuum light curves for each 
AGN, including the low-variability objects shown in Figure 

M 

Table |2] also presents measurements of f\( 5100 A), the 
AGN continuum flux density at 5100 A in the AGN rest 
frame. This quantity is most often used to determine AGN lu- 
min osities for use in the radius-luminosity relationship (e.g., 
iBentz et akll2009al) . This wavelength corresponds to a local 
minimum in Fe II flux, making it a good choice for measure¬ 
ment of the featureless continuum level in total-flux spectra. 
We used the AGN power-law component of the fit to the mean 
spectrum to measure f\ (5100 A) for each AGN. The formal 
fitting uncertainty on the power-law component flux for the 
mean spectrum is typically very small (a few percent) and not 
a good measure of the actual uncertainty. From examination 
of the fitting results using all three Fe II templates, we find 
that the rms scatter among the power-law component fluxes at 
5100 A for the three fits ranges from 4% to 13% of the mean 
value. For uniformity, we adopt an estimated 10% uncertainty 


Mrk 40 



HJD - 2450000 


Fig. 7.— Light curves of Mrk 40. The y-axis units for this and subsequent 
light curve plots are 10 15 erg cirT 2 s -1 A 1 for the [4-band continuum flux 
density and 10 12 erg cm ~ 2 s 1 for the emission-line fluxes. Plotted error bars 
represent the propagated uncertainties from the spectral extractions only, as 
listed in Table |4] and do not incorporate the estimated residual flux normal¬ 
ization scatter as determined from the [O III] emission line. 

on the 5100 A featureless continuum flux at this wavelength 
resulting from differences among the Fe II templates. We 
also combine with this (by addition in quadrature) an esti¬ 
mate of the overall photometric uncertainty on the continuum 
flux density given by <5([0 III])//([O III])x f\( 5100 A), where 
the [O III] flux and uncertainty values are listed in Table [2] 
Improved values of /\(5100 A) can be determined using HST 
V-band images to remove the host-galaxy contribution more 
accurately from the total spectroscopic flux densitvPl 

For the strongly variable AGNs, we generally find higher 
values of /\ al for the higher-order Balmer lines: F var (H<5) > 
F var (H 7 ) > F var (H/3) > F val -(Ha:). This overall trend is consis¬ 
tent with results from previous work fe.g.. IBentz et al.l2010al) . 
and is i n acco rd with e xpecta tions from photoionization mod¬ 
eling dKorista & Goadll2004l) . The He II A4686 line is by far 
the most strongly variable line in the optic al spectr um, follow¬ 
ing the trend first seen in NGC 5548 by iPeterson & Ferlandl 
(1986). For the six objects having He II light curves, the ratio 
/\ a ,(He II)/F var (H/3) has a mean value of 3.7. The strong re- 
sponsivity and prompt response of He II is easily visible in the 
light curves: for example, in Mrk 40, the He II emission drops 
very steeply following the decline in continuum luminosity at 
about HJD 5685, while the decline in the Balmer-line fluxes 
is much shallower and more gradual (Figure IT}. 

32 New images of Mrk 50. Mrk 704, Mrk 1511, and Zw 229-015 will be 
obtained with WFC3 in the F547M filter as part of HST Cycle 22 program 
GO-13816 (PI: Bentz) along with 10 other recently reverberation-mapped 
AGNs compiled from other programs. 
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Fig. 8 .— Light curves of Mrk 50, including observations taken during FIG. 10.— Light curves of Mrk 1511. 

January-March before the start of the main campaign. 
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FIG. 9.— Light curves of Mrk 279. 


FIG. 11.— Light curves of Mrk 704. 



































17 


NGC 4593 
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FIG. 12.— Light curves of NGC 4593. 


Zw 229-015 
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Fig. 13.— Light curves of Zw 229-015, including the three data points 
taken after the end of the main campaign. 


While it is possible to measure reverberation lags based on 
the data presented here for the highly variable AGNs in our 
sample, we find that for some objects the cross-correlation 
measurements show evidence for an excess and spurious sig¬ 
nal at zero lag, caused by correlated flux-calibration errors in 
the continuum and emission-line light curves. Thus, we defer 
the measurement of lags to a later paper in this series, in which 
our emission-line light curves will be combined with V -band 
light curves which will be presented separately. The V -band 
light curves generally have a higher cadence and higher S/N 
than the U s light curves. Determination of virial masses for 
the black holes in these AGNs will be done using lags mea¬ 
sured against the V-band light curves. 

6. BROAD H/3 WIDTHS AND VELOCITY SHIFTS 
6.1. Instrumental Resolution 

To quantify the broad-line widths, it is necessary to cor¬ 
rect measurements for the instrumental broadening in obser¬ 
vations with the 4"-wide slit. The spectral resolution for 
observations of the intrinsically compact BLR as observed 
through a 4" slit will be seeing-dependent, and will vary 
slightly from night to night due to variations in spectrograph 
focus. However, the impact of instrumental broadening on 
the broad emission lines in these AGNs is very small, and we 
use a single representative value for the instrumental broaden¬ 
ing for simplicit y. For the b l ue-side data we follow the same 
method used by IBentz et ahl (l2009bl) and iBarth et afl (1201 lal) 
to estimate the instrumental FHWM by using observations of 
an AGN and an arc lamp in 2" and 4" slit widths. Mrk 40 
was observed on the night of 2011 May 4 in both of these slit 
widths. For the spectral region around H (3, we use the Cd I 
A5085 A line as the nearest reference calibration line. From 
the 2" observations of the AGN and the comparison lamp, 
we subtract the width of the Cd I line in quadrature from the 
observed [O III] width to derive an estimate of the intrinsic 
[O III] line width of FWHMw 4.9 A. Then, the instrumen¬ 
tal broadening for the 4" slit is taken to be the difference in 
quadrature between the [O III] line width as observed through 
the 4" slit (FWHM = 8.1 A) and the intrinsic [O III] width. 
This gives an instrumental FWHM of 6.45 A or 380 km s -1 , 
and an instrumental dispersion of cri nst m 162 km s -1 . 

For the red-side data, there is no strong narrow line that can 
be used as a reference for comparison between the 2" and 4" 
observations, and instead we use the widths of comparison 
lamp lines as observed through the 4" slit as a rough esti¬ 
mate of the instrumental broadening, obtaining an instrumen¬ 
tal FWHM of 7.7 A or ~ 330 km s -1 for the region around the 
Ha line. These instrumental widths are much smaller than the 
broad emission-line widths in the AGNs, and the corrections 
to the observed line widths are very small. 

6 .2. Broad 118 Widths in Mean and RMS Spectra 

The width of the broad H/3 line is a primary quantity used 
for determination of BH masses from reverberation data, and 
it is of interest to measure both the FWHM and the disper¬ 
sion (or second moment) of the H/3 line profile in both the 
mean and the rms spectra when possible. We use the results 
of our spectral decompositions to measure these width param¬ 
eters for the broad H/3 line after isolating it from the con¬ 
tributions of other emission-line and continuum components. 
For each nightly spectrum, we subtract all model components 
other than broad H/3 from the data, leaving the broad H/3 pro¬ 
file as the residual. From this time series of residuals, we 
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Light curves of the remaining objects in the sample. In each panel, the upper and lower subpanels show the U s and H/3 light curves, respectively. 


FIG. 14.— 

construct mean and rms spectra, and measure widths empiri¬ 
cally from the mean and rms line profiles as described below. 
The mean spectrum has a continuum level of zero (modulated 
by residual noise from the spectral decompositions), while the 
continuum level in the rms spectrum is nonzero due to photon¬ 
counting noise and residual errors in the nightly spectral de¬ 
compositions. This noise continuum level in the rms spectra is 
generally flat and featureless outside the H/3 profile (as seen in 
the bottom panel of Figure [6| and can be subtracted off with 
a simple linear fit to regions on either side of the line. For 
each AGN we selected broad continuum windows redward 
and blueward of the H/3 line to fit and subtracted this contin¬ 
uum from the rms spectrum, avoiding any residual spikes at 
the wavelengths of [O III] or other strong features. 

Line width parameters were measured following methods 
similar to those of lPeterson et al.l (120041) . The profile FWHM 
was measured by starting at the profile peak and moving red- 
ward until reaching a flux equal to half the peak value (using 
linear interpolation between pixels on the profile shoulder to 
determine the wavelength at half maximum), and similarly 


moving blueward until reaching half of the peak flux, and de¬ 
termining the separation in Angstroms between the two half¬ 
peak points. To measure the line dispersion <7ii ne for either 
the mean or rms profile, we first determine the line centroid 
wavelength Ao: 

\ E ^ifi /ox 

0 Eft ’ 

where f, is the flux density at pixel i and the summation is car¬ 
ried out over all pixels in the profile. Then, the line dispersion 
is calculated according to 

2 _ ( EEfi \ _ x2 / 4 x 

lme V Ef ) °' 

The endpoints of the summation region are determined visu¬ 
ally for each AGN by identifying the location where the line 
profile meets the surrounding continuum level. 

In order to assess the uncertainty in the measured width pa¬ 
rameters, we carry o ut a Monte Carlo procedu re, following 
iPeterson et al.l (120041) and iBentz et afl (l2009bl) . by creating 
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10 4 modified realizations of the mean and rms profiles and 
measuring the width parameters for each one. Each realiza¬ 
tion is constructed from the time series of spectra by randomly 
selecting N individual nightly spectra (where N is the total 
number of spectra in the time series), allowing duplications, 
calculating the mean and rms from this revised set of spec¬ 
tra, and measuring the broad H/3 widths from each. For mea¬ 
surements of the rms profile, we allow the endpoints of the 
continuum fitting region in each realization to vary randomly 
within a range of ±10 A from our initial selections, and for 
measurements of line dispersions from either the mean or rms 
profiles, we allow the endpoints of the summation region in 
each realization to vary randomly within ±10 A from our ini¬ 
tial selections. The final line width and uncertainty are taken 
to be the mean and standard deviation of the distribution of 
measurements over these 10 4 iterations. 

Finally, the line widths are converted to velocity units, and 
the instrumental FWHM or a are subtracted in quadrature 
from the measured FWHM or a values (respectively) to cor¬ 
rect for instrumental broadening. The results of these mea¬ 
surements are reported in Table [3 Width parameters for 
the rms profile are only reported for the five AGNs having 
fvarf H/i) >0.1 and for Mrk 279. For the remainder of the 
sample, the rms spectra contain very little flux as a result of 
the overall low variability amplitude, and the Monte Carlo er¬ 
ror analysis procedure produces many iterations in which the 
FWHM or (Tiine are not well-defined. Widths measured from 
the rms profiles for these low-variability AGNs have large un¬ 
certainties and would not be particularly useful for deriving 
estimates of black hole masses, despite the fact that some of 
these low-variability AGNs appear to show clear signal at H/3 
in the rms spectra (Figured- 

An important caveat for the line width measurements, par¬ 
ticularly for the line dispersions, is that the red wing of H/3 
is often blended with a substantial contribution of Fe II and 
possibly He I as well. As described above, the exact amount 
of this contamination is difficult to determine, and the end¬ 
point of the red wing of the line is usually not well deter¬ 
mined. Since the line dispersion values are highly sensitive 
to the amount of flux in the line wings, broad H/3 line disper¬ 
sions are particularly susceptible to systematic errors resulting 
from this deblending. We note that in most past reverberation¬ 
mapping work, no deblending was attempted at all, potentially 
leading to overestimates of H/3 line dispersions. The FWHM 
may be modestly affected by this red-wing contamination as 
well. 

6.3. Broad H/3 Width and Velocity Centroid Variations 

Width parameters (FWHM and erii ne ) were also measured 
from each nightly spectrum for the high-variability objects in 
order to search for time-dependent changes in profile width, 
following the same methods described above. These widths 
were measured from the spectra of the broad H/3 compo¬ 
nent alone, constructed by subtracting all other model com¬ 
ponents (including the narrow H/3 component) from each to¬ 
tal flux spectrum. To assess the uncertainty in FWHM and 
Cline, we created 1000 modified realizations of each spectrum 
by adding Gaussian random noise with the noise dispersion 
equal to the amplitude of the error spectrum at each wave¬ 
length, and then measuring the FWHM and crij ne from each 
modified spectrum. The uncertainties on FWHM and erij ne 
were taken to be the standard deviations of the distributions 
of values measured from this procedure. In general, the line 


Mrk 40 
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Fig. 15.— Broad H/3 parameters as a function of time for Mrk 40. The 
top panel shows the same H/3 light curve as in Figure [7] in units of 10 1 5 
erg crrT 2 s . Middle panels show the broad H/3 FWHM and line dispersion 
£T]j ne , and the bottom panel shows the velocity shift AV measured relative to 
the mean spectrum. 

widths show an anticorrelation with flux; §(7] discusses this 
trend further. 

Using the same spectra of the broad H/3 component, we 
measured the broad-line velocity centroid shift relative to 
the broad H/3 component of the mean spectrum, using a 
cro ss-correlati o n met hod based on the tec hniques described 
by IShen et akl (12013l) and lJu et al.l (120131) in the context of 
searching for binary black hole orbital motion in quasars. The 
wavelength centroid of the broad H/3 component is a param¬ 
eter determined by our spectral fitting procedure, but the fit¬ 
ted line centroid can be strongly influenced by asymmetric 
wings on the line profile, particularly in the H/3 red shelf. 
In contrast , the c ross- correlation line shift as measured by 
IShen et al.l dMl an d[Ju et al.l (120131) is a better measure of 
an overall wavelength shift of the line as a whole, and is less 
sensitive to noise in the high-velocity wings. The velocity 
shifts are discussed in §j8] 

The nightly width and velocity shift parameters are dis¬ 
played in Figures fBUTOl for Mrk 40, Mrk 50, Mrk 279, 
Mrk 1511, and NGC 4593, all of which show distinct time- 
dependent trends in line profile width and/or velocity. For the 
remainder of the sample, no clear signal was detected beyond 
random scatter in the measurements. 

7. RELATIONSHIPS BETWEEN BROAD-LINE WIDTH AND 
LUMINOSITY 

Figures 1201423 display the relationship between broad H/3 
width (FWHM or <7ij ne ) and broad H/3 flux. The narrow H/3 
component flux from the mean spectrum decomposition has 
been subtracted from the total H/3 light curve in each case. 
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TABLE 5 

Broad H/3 Component Widths in Mean and RMS Spectra 



Mean Spectrum 

RMS Spectrum 

Galaxy 

FWHM 

^line 

FWHM 

^line 


(km s -1 ) 

(km s' 1 ) 

(km s -1 ) 

(km s _1 ) 


Mrk 40 
Mrk 50 
Mrk 141 
Mrk 279 
Mrk 486 
Mrk 493 
Mrk 504 
Mrk 704 
Mrk 817 
Mrk 841 
Mrk 1392 
Mrk 1511 
NGC 4593 
PG 1310-108 
Zw 229-015 


2021 ± 17 
4101 ±56 
5129 ±45 
4099 ±43 
2106±8 
1717 ±25 
1927 ± 7 
6664 ±66 
5593 ±32 
6065 ±34 
5421 ±31 
4154 ±28 
4264 ±41 
3422 ±21 
3705 ±203 


1058 ±9 
2024 ± 31 
2280 ±21 
1821 ±13 
1332 ±14 
1128 ±41 
1058 ±14 
2911 ±61 
2718 ± 45 
2794 ±48 
2103 ±17 
1828 ± 12 
1925 ±38 
1823 ±20 
1747 ±56 


1688 ±26 
3355 ±128 

3306 ±338 


3236 ±65 
3597 ±72 

1789 ±93 


740 ±17 
2020 ± 103 

1778 ±71 


1506 ±42 
1601 ±40 

1609 ± 109 


NOTE. — Widths are measured from the mean and rms of the broad 
H/3 component spectra, as described in the text. 
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Fig. 16.— Broad H (3 parameters for Mrk 50. Only the time span of the 
main campaign is shown, because the data points from earlier observations 
are very noisy. 

to obtain a light curve of the broad H/3 line alone. The flux 
uncertainties have been expanded by including the estimated 
uncertainty in spectral scaling determined from the [O III] ex¬ 
cess scatter, using the results listed in Table [2] In the case of 
Mrk 50 (T'igure fTil we include only data points from the main 
observing campaign, since the observations taken prior to the 
start of the campaign are much noisier. Mrk 279 is excluded 
since the broad H/3 width could not be measured accurately, 
as described previously. 


Fig. 17.— Broad H/3 parameters for Mrk 279. The nearly constant H/3 
width for this AGN is an artifact of the constraint that had to be imposed in 
order to prevent the line width from reaching unphysically large values in the 
fits. 

Although there is substantial scatter in some of the plots, in 
each case we find an overall anticorrelation between the broad 
H/3 width and flux. To quantify this relationship, we use the 
Bayesian linmix_err regression method of lKellvl (120071) . 
fitting a model of the form log[AV(H/3)] = /31og[/(H/3)] + a, 
where AT represents either FWHM or <7ij ne . The slopes de¬ 
termined by the linmix_err regression are shown in each 
figure and range from /3 = -0.07 to -0.28. 

Similar trends between broad H/3 width and luminosity 
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8. — Broad H/3 parameters for Mrk 1511. 
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FIG. 19.— Broad H/3 parameters for NGC 4593. 


have been seen previously in high-quality rever beratio n map¬ 
ping d atasets. Examples inc lude NGC 5548 dDennev et alJ 
l2009al) and Mrk 40 (Arp 151: lPark et alj|20l2ah . The anticor¬ 
relation between broad H/3 width and luminosity arises as a 
result of the relationship between the ionizing flux and the lo¬ 
cal responsivity (or reprocessing efficiency) of BLR gas: the 
H p responsivity is highest in the outer port ions of the B LR 
where incident continuum fl uxes are lower dKorista & Goadl 
l2004t iGoad & Koristall20l4h . An increase in ionizing con¬ 
tinuum luminosity will therefore lead to an increase in the 
emissivity-weighted BLR radius; this effect is known as 
“breathing”. For a BLR having a radial velocity gradient with 
the highest-velocity gas at small radii close to the black hole, 
breathing will preferentially boost the flux of the low-velocity 
line core relative to the high-velocity wings of the broad line, 
making the line profile narrower. 

BLR breathing has usually been considered in the context 
of large luminosity changes over sufficiently long timescales 
that emission-line lags can be seen to vary betwe en widely 
separated epochs of monitoring (e.g.. lO’Brien et al.l 11995 : 
Korista & Goadl 120041 ICackett & Hornd l2006t IBentz et alJ 
20071) . The H B width-lumin osity anticorrelati on seen in our 
data (and in recent work by iPark et al.l:2012 a) further illus¬ 
trates that BLR breathing occurs even on the shortest observed 
timescales (days to weeks) in response to modest variations 
in the AGN continuum. While our data do not directly reveal 
short-timescale changes in lag or BLR radius, the relationship 
between H/3 width and luminosity implies that the BLR ra¬ 
dius responds very rapidly to continuum flux changes. A fur¬ 
ther consequence of this finding is that the reverberation trans¬ 
fer function will be nonstationary even on timescales shorter 
than a single typical monitoring campaign. Measu red trans¬ 
fer functions (IBentz et al.ll2010bc iGrier et afll2013al) will thus 
represent temporally smeared averages of intrinsically time- 
varying functions. 

Since broad-line flux variations are driven by continuum 
variations, it is also of interest to examine the relationship 
between broad-line width and continuum flux. However, the 
time delay between continuum and line variations introduces 
substantial additional scatter to the relationship, and in our 
data, the broad H (3 widths show a much clearer relation¬ 
ship with H/3 flux than with continuum flux. The temporal 
mismatch between line width and continuum flux measure¬ 
ments can be partially alleviated by selecting continuum and 
emission-line measurements that are sep arated by the mea¬ 
sured emission-line lag time (IPark et al.ll2012ali . Unfortu¬ 
nately, the frequent weather-related gaps in our light curves do 
not leave enough matched pairs of data points to obtain clear 
results. Another possibility would be to interpolate the data 
to obtain matched pairs of line widths and continuum fluxes 
after correcting for the broad-line lag time, but we will defer 
that investigation in order to use the higher-cadence V -band 
photometric light curves instead of the spectroscopic contin¬ 
uum measurements presented here. 

The behavior of broad-line widths in response to luminos¬ 
ity variations has important ramifications for the accuracy 
of single-epoch spectroscopic determinations of black hole 
masses. These are derived by using the AGN continuum 
luminosity as a proxy for the BLR radius and computing a 
virial product based on VP = L 7 ( AV) 2 /G, where 7 represents 
the slope of the radius-luminosity relationship. Recent work 
finds rBLR oc L ^ 5 ^ 0 ' 04 for the full sample of reverberation- 
mapped AGNs, where L 5100 is the AGN continuum luminosity 
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Fig. 20.— Relationship between broad H ft width and flux for Mrk 40. In 
each panel, the solid line shows the best fit to the relationship with slope ft. 
Line widths (FWHM and <rij ne ) are in km s -1 , and flux units are 10 -15 erg 

-2 -1 

cm s . 

at A = 5100 A dBentz et al.l2013l) . The luminosity of the broad 
Balmer emission lines can also be used as a proxy for contin¬ 
uum luminosity dGreene & Hdl2005 i. The accuracy of these 
single-epoch methods for mass estimation relies implicitly on 
the assumption that line width and luminosity must be anticor¬ 
related in a given AGN, such that the derived masses will not 
change dramatically as the luminosity varies. Our results pro¬ 
vide add itional reassur i ng evid enc e that this is gener ally the 
case; see lDennev et al.1 (I2009al) and lPark et alJ (120 12 all for de¬ 
tailed discussion of the impact of variability on single-epoch 
mass determinations. If ielr oc L ( ^ f 00 and the H/3 luminos¬ 
ity scales linearly with continuum luminosity, then an AGN 
maintaining a constant virial product would exhibit a line 
width-luminosity relationship having slope /3 = -0.25, similar 
to the slopes observed in Mrk 40, Mrk 1511, and NGC 4593. 
However, it is important to note that this naive prediction for 
the slope of the H/3 width-luminosity relationship is an over¬ 
simplification; for example, observations of NGC 5548 show 
that it exhibits more complex behavior. The broad H/3 flux in 
NGC 5548 does not scale linearly with continuum flux (the 
“intrinsic Baldwin effect”; [Gilbert & Petersonll2003l) . and the 
relationship between BLR radius and continuum luminosity 
measured from long-term monitoring of NGC 5548 exhibits 
a steeper slope than the radius-luminosity relationship for the 
reverberation-mapped sample as a whole ( iKilerci Eser et al.l 
l20l5l) . 

Recently. IGuo & Gul (120141) carried out line width measure¬ 
ments for a sample of 60 quasars having at least six epochs 
of SDSS spectroscopy, and they examined the correlation be¬ 
tween line FWHM and luminosity for C IV, Mg II, H/3, and 
Ha. Interestingly, 75% of their sources show a positive cor¬ 
relation between broad-line width and luminosity, and only 
three sources exhibited an anticorrelation at the > 90% confi¬ 
dence level. Another example of a positive width-luminosity 
correlation for the Balmer lines was recently found in PG 
1613+658 (lZhangll2014l) . using data from the Wise Observa¬ 
tory quasar monitoring campaign dKasni et al.l 12000 1. This 
behavior is very different from that seen in our sample, and 
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Fig. 21.— Relationship between broad H/3 width and flux for Mrk 50. Only 
data points from the main campaign are included. 
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Fig. 22.— Relationship between broad H/3 width and flux for Mrk 1511. 
is inconsistent with expectations from the scenario of BLR 
breathing determined by radial gradients in responsivity (with 
higher responsivity at larger radius) as described above. It 
will be interesting to examine this trend for larger samples of 
quasars over long time baselines as more multi-epoch spectra 
become available for quasars over longer time baselines. If 
broad-line widths in high-luminosity quasars respond to con¬ 
tinuum variations in a manner strongly different from what we 
observe in our sample of Seyfert galaxies, this would poten¬ 
tially compromise the reliability of single-epoch black hole 
mass estimates derived from quasar spectra. 

8 . IMPLICATIONS OF BROAD-LINE VELOCITY SHIFTS FOR 
BINARY BLACK HOLE SEARCHES 

The spectral decompositions make it possible to measure 
night-to-night changes in the broad H/3 velocity centroids. 
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NGC 4593 



Fig. 23.— Relationship between broad H/3 width and flux for NGC 4593. 


Our results for Mrk 40, Mrk 50, Mrk 279, Mrk 1511, and 
NGC 4593 illustrate that systematic shifts of ~ 100 km s _1 , 
correlated with flux changes, are common and occur over 
timescales similar to the continuum variability timescale. The 
largest shifts are seen for NGC 4593, in which the broad H/3 
centroid changed by 266 ±11 km s -1 over a 30-day interval, 
as seen in Figure [19] 

Broad-line velocity shifts are expected in any AGN for 
which the reverberation transfer function is asymmetric 
about the line center, which will occ ur when radial mo- 
tions are present in the BLR (e.g.JBlandford & McKeelll982l: 
Canriotti et al.l[l982j: [Antonucci & Coheni 1983: Ulrich et all 


1984tlStirne etdll988tlRosenblatt & Malkanll 19901) . The re¬ 


verberation process can be described in terms of a convolu¬ 
tion of the continuum light curve C(t) by a velocity-dependent 
transfer function. 


L(v r ,t)= / *(v n T)C(t-T)dT, (5) 

Jo 

where the transfer function \E'(vy,t) is a function of line- 
of-sight velocity v r and of time delay t, and L(v r .t) 
gives the resulting time- dependent emission-line profile 
( Blandford & McKedlT98l PeTei>on IMlIl ). If BLR clouds 
are in circular orbits around the black hole, as in a circular 
Keplerian disk, then the transfer function will show short lags 
at high velocities on both the blueshifted and redshifted sides 
of the line profile, and at small line-of-sight velocities there 
will be a broad distribution of lag time. If the BLR gas is 
dominated by radial motions, however, the transfer function 
will be strongly asymmetric about the line center. For ex¬ 
ample, in a BLR dominated by radial inflow, the redshifted 
gas on the near side of the BLR will have near-zero lag, 
while the blueshifted gas on the far side of the BLR will have 
the longest lags. Representative tra nsfer functions for sim ¬ 
ple kinemat i c mode l s are shown bvlWelsh & Home! (119911). 
IBentz et all (1 201 Obi) . iGrier et aD (1201 3al) . and IPancoast et all 
(I2014bl) . 

In addition to the imprint of BLR kinematics on the trans¬ 
fer function, the transverse Doppler shift and gravitational 


redshift in the BLR will produce redward asymmetries in 
thejransfer function, particularly in the high-velocity wings 
dGoad et afli2012l) . The combination of transverse Doppler 
shift and gravitational redshift, together with any radial mo¬ 
tions, should ensure that all transfer functions are at least 
modestly asymmetric in structure even if the BLR cloud mo¬ 
tions are predominantly rotational. Absorption or scattering 
within the BLR may also alter the shape of 'IT 17, r) in an 
asymmetric fashion. 

With the best-quality reverberation-mapping data obtained 
over the past several years, it has become possible to ex¬ 
amine the velocity-dependent structure of the transfer func¬ 
tion in several low-redshift AGNs by measuring the la g as a 
function of veloc i ty across the fine profile (lKo llatschn vll2003l : 

IBentz et alJ2008tlDennev et al.l2010llBarth et all201 lbl) . and 

by reconstru cting the full two-dimensional shape of the trans¬ 
fer function (IBentz et al.ll2010bl : lGrier et al.li2013 a). Dynam¬ 
ical modeling of reverberation data has also begun to yield 
constrai nts o n the transfer function shape (iPancoast et al.l 
I 20 T 2 L l20T4bh . A primary result that has emerged from these 
investigations is that the H f3 transfer functions for nearby 
Seyfert galaxies exhibit a broad diversity of structures: some 
objects have nearly symmetric transfer functions consistent 
with expectations for a Keplerian disk, while others have 
strong asymmetries suggestive of inflowing or outflowing mo¬ 
tions. 

As an illustration of why an asymmetric transfer function 
will lead to broad-line velocity centroid changes, consider 
a radially inflowing BLR. If the continuum luminosity in¬ 
creases, the prompt red-wing response will shift the line cen¬ 
troid redward initially until the blue wing response from gas 
on the far side of the BLR is seen at a later date. Similarly, a 
decrease in continuum luminosity will initially shift the line 
centroid blueward. The opposite behavior would be seen in 
an outflowing BLR. Line centroid variability will depend on 
the continuum light curve shape and variability amplitude, the 
broad-line width, and the degree of asymmetry in the transfer 
function. For a specified transfer function, the direction of 
line profile shifts will depend on the recent gradient in the 
continuum light curve rather than on the absolute luminosity 
of the continuum. Either blueshifts or redshifts can occur at 
a given continuum flux level, depending on the recent history 
of continuum variations. 

These considerations are relevant for binary black hole 
searches because one proposed method for identification of 
binary black holes involves searching for temporal variations 
in th e line-of-sight velocity of broad em ission lines in quasars. 
(See iPopovid 120 1 2l iDotti et aD 120121 for recent reviews of 
binary black hole phenomenology, optical spectroscopic sig¬ 
natures, and search techniques). If one member of a binary 
black hole is active and possesses a broad-line region, then 
the orbital motion of the binary will cause the broad emission 
lines to oscillate in velocity over an orbital period. A sub¬ 
stantial amount of effort has been invested recently in search¬ 
ing for broad-line velocity changes in quasars having two or 
more epochs of spectroscopic observations. One strategy is to 
search for objects having very large broad-line velocity off¬ 
sets (i.e., |AV| > 1000 km s -1 relative to the host galaxy), 
and then obtain follow-up spectroscopy to check whether the 
broad-line velocity drifts over time in a m anner con siste nt 
with expectations f or orbital motion (lEracleous et al. 1 201 2t 
iDecarli et all 120131 : iTsalmantza et al.l 1201 lb iLiu et al.l 20141) . 
An alternate method is to start with the general quasar popula- 
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tion, rather than restricting to objects showing strong velocity 
offsets initially, and to search for changes in broad-line cen¬ 
troid velocity between two widely separated epochs of spec¬ 
troscopic observations (IShen et alJl2013tlJu et al.ll2013l) . The 
discussion in this section is primarily (but not exclusively) rel¬ 
evant to this latter method. 

Such searches are still in their initial stag es. In recent 
searches using SDSS data, IShen et aTI (120131) examined H/3 
velocity shifts in a sample of 1347 quasars having two or more 
epochs of spectroscopy, finding 28 objects with significant 
velocity shifts (of typically ~ 100-300 km s -1 ) and identi¬ 
fying se ven of these as being promising candidate binaries. 
iJu et ah . (12013t) carried out a similar search using Mg II, find¬ 
ing seven candidate sub-pc binaries with velocity changes of 
> 280 km s -1 from a starting sample of 1523 quasars having 
spectroscopy with S/N > 10 at two epochs. For these candi¬ 
dates, continuing long-duration follow-up spectroscopy will 
be required in order to rule out alternative explanations and 
conclusively demonstrate binary motion. 

In discussion of these search efforts, it is generally recog¬ 
nized that broad-line velocity changes can happen as a result 
of structural changes in the BLR occurring over_a dynami¬ 
cal timescale (years to decades) or longer (IShen et al.ll2013l) . 
and such changes have been seen in multi-yea r mo nitoring ob- 
serva ti ons of NGC 5548, for example (Wanders & Peterson! 
Il996t iSergeev et alJ 120071) . The appearance or fading of 
hotspots in the BLR has also been suggested as a possi¬ 
ble mechanism to induce line profile changes and apparent 
velocity shifts (iDecarli et al.l 120131) . However, asymmetric 
reverberation seems to have been overlooked as a possible 
cause of the detected velocity shifts in these binary black 
hole searches. Our data provide a clear demonstration that 
broad-line velocity centroids can vary significantly on short 
timescales solely as a result of reverberation, and do not pro¬ 
vide direct evidence for bulk acceleration of the BLR. 

Based on our small sample of nearby, low- to moderate- 
luminosity AGNs, it is difficult to extrapolate our results to 
make any useful predictions for the distribution of broad-line 
velocity shifts expected as a function of timescale for quasar 
surveys. We have carried out simple simulations using asym- 
metric transfer func tions similar to those observed in Arp 151 
(Bentz et al.l2010bt) or NGC 5548 (Pancoast et al.ll2014bl) and 
find that velocity shifts of > 300-400 km s _I can occasion¬ 
ally occur for plausible AGN variability amplitudes. How¬ 
ever, the range of transfer function shapes and asymmetries is 
completely unknown for the quasar population. For a given 
AGN variability model and transfer function model, it would 
be straightforward to predict the expected distribution of ve¬ 
locity changes as a function of timescale, but such an effort is 
beyond the scope of this paper. 

Here, our goal is simply to call attention to the fact that re¬ 
verberation does indeed produce significant velocity shifts on 
short timescales as a result of the BLR response to continuum 
fluctuations. This has been understood since the earliest re¬ 
verberation mapping studies, and our spectral decompositions 
provide a way to examine broad H/3 velocity shifts in Seyfert 
1 nuclei with high-cadence sampling over a 2.5-month base¬ 
line. Considering that we find one object showing a velocity 
change of > 250 km s -1 in our small sample over this brief 
monitoring duration, it seems entirely plausible that, in large 
AGN samples observed over long time baselines, the tail of 
the distribution of reverberation-induced velocity shifts could 
extend to substantially larger velocities than the values we ob¬ 


serve here. It is difficult to exclude the possibility that most or 
all of the recently identified candidate binaries might simply 
be exhibiting velocity changes due to this relatively prosaic ef¬ 
fect, rather than exhibiting acceleration due to binary orbital 
motion. In large samples of quasars with two-epoch spec¬ 
troscopy, selection of targets with the largest velocity shifts 
will preferentially select objects from the tail of the distri¬ 
bution of large reverberation-induced velocity shifts, and we 
suggest that the majority of candidate binaries selected via 
this method will most likely be false positives. 

This reverberation effect will also contribute random ve¬ 
locity noise or jitter that could impede the detection of gen¬ 
uine binaries (particularly wide-separation or nearly face-on 
binaries), rather analogous to the impact of stellar activity on 
the detectability of exoplanets via the radial-velocity method. 
Furthermore, the radial-velocity curve for NGC 4593 (Fig¬ 
ure [l9l> illustrates another potential concern. In this case, the 
velocity curve could be mistaken for a near-complete cycle 
of sinusoidal oscillation, when in reality the line centroid is 
simply responding to a down-and-up variation in the contin¬ 
uum luminosity. In order to convincingly demonstrate binary 
orbital motion and rule out reverberating impostors similar to 
this example, it will be necessary to monitor broad-line veloc¬ 
ity variations over at least two orbital cycles, with sufficient 
cadence to sample the orbit well, and also to demonstrate that 
the velocity variations are uncorrelated with continuum flux 
variations. The expected orbital periods for quasars hosting 
su permassi ve black hole binaries span years to decades (e.g., 
IShen et aUl2013l) . thus confirmation of black hole binaries by 
this method will require a very long-term effort. 

9. CONCLUSIONS 

We present spectroscopic data, light curves, and other mea¬ 
surements from a 69-night reverberation mapping campaign 
carried out at the Lick 3 m telescope in Spring 2011 for a sam¬ 
ple of 15 low-redshift Seyfert 1 galaxies. Although our ob¬ 
serving run was hampered by worse-than-average weather, we 
obtained data of excellent quality for a subset of our sample, 
suitable for measurement of light curves for multiple emission 
lines. Our primary results can be summarized as follows. 

1. Expanding on our team’s recent work, we have car¬ 
ried out multicomponent decompositions of the nightly 
blue-side spectra from our campaign. The decomposi¬ 
tions enable measurement of light curves for blended 
features including Fe II and He II, and measurements of 
H/3 light curves and widths are improved by deblending 
H/3 from other contaminating emission-line and contin¬ 
uum components. Obtaining a time series of H/3 line 
profiles without contamination from other overlapping 
features is also important for measurement of velocity- 
resolved reverberation signals and transfer functions, 
and for application of dynamical modeling methods to 
derive black hole masses (Pancoast et alJl2014bl) . In 
Appendix O below, we demonstrate that fitting and re¬ 
moval of continuum components prior to constructing 
the rms spectrum is also necessary in order to avoid a 
potentially severe bias that can occur when measuring 
line widths from the rms spectrum. Spectral decompo¬ 
sitions should be carried out as a routine part of rever¬ 
beration mapping data analysis, provided that the S/N 
of the data is sufficient for multicomponent fits to yield 
useful results. 

2. General variability properties for our sample are consis- 
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tent with trends seen in previous reverberation-mapping 
programs, including higher responsivity for the higher- 
order Balmer lines, and very high responsivity for He II. 
While cross-correlation lags can be measured between 
the emission-line and continuum light curves presented 
here, lags measured against V-band photometric light 
curves will have better reliability as a result of the 
higher S/N of the photometric data and the absence 
of correlated flux-calibration errors between the pho¬ 
tometric and spectroscopic data. 

3. The anticorrelation between broad-line width and flux 
observed in four objects demonstrates that the BLR 
“breathes” (i.e., the emissivity-weighted BLR radius 
increases as a function of AGN luminosity) even on 
short timescales of days to weeks in response to contin¬ 
uum variations. This implies that reverberation transfer 
functions will be nonstationary even during the course 
of a single monitoring campaign. This breathing effect 
is responsible for the relationship between broad-line 
profile widths measured in mean and rms spectra, as 
we illustrate by means of simulations in Appendix [B] 

4. In NGC 4593 we observe a broad H/3 velocity shift 
of ~ 250 km s -1 over the course of one month, and 
other AGNs in our sample exhibit H/3 velocity shifts 
of order 100 km s -1 over similar timescales. Such ve¬ 
locity shifts caused by asymmetric reverberation repre¬ 
sent an important source of confusion noise for binary 
black hole searches based on multi-epoch quasar spec¬ 
troscopy. When only two (or a few) observations of 
a given AGN are available, line-profile centroid shifts 
due to velocity-dependent reverberation can mimic the 
expected acceleration signature of binary orbital mo¬ 
tion, producing a spurious identification of a candidate 
binary. Typical velocity changes caused by asymmet¬ 
ric reverberation will be up to a few hundred km s -1 , 
similar to the velocity shifts detected in many of the 
recently identified candidate binaries (IShen et al.ll2013l: 
lJu et al.ll2013h . In general, reverberation will produce a 
velocity-jitter noise floor that could limit the detectabil¬ 
ity of binary black holes in spectroscopic searches for 
radial-velocity oscillations. 

Additional papers in this series will present V-band pho¬ 
tometric light curves, measurements of broad-line lags and 
velocity-resolved reverberation using both cross-correlation 
and JAVELIN ( Zu et al.l 12011 1 methodology, black hole 
virial mass estimates, and BLR d ynamic al modeling us- 
ing the methods des cribed by iPancoast et al.l (1201 111 and 
iPancoast et al.l (120 Hall . Following the completion of this pro¬ 
gram, our reduced and scaled spectra will be made available 
to the community for archival use FI 
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Fig. 24.— Kast spectra of objects observed at the start of the campaign but not selected for continued monitoring. 


APPENDIX 

A. ADDITIONAL SPECTRA 

During the first few clear nights of our Lick program (2011 March 29-31), we observed several additional galaxies that were 
initially selected as possible reverberation targets based on SDSS or other spectra or on AGN classifications from the literature, 
but which were subsequently discarded from our monitoring sample. Some of these objects were too faint for regular monitoring, 
and others appeared to be in a low flux state compared with previous observations. In one case we did not detect any emission 
lines in a galaxy previously classified as a Seyfert 1. Interestingly, three AGNs selected from SDSS based on their strong broad 
H [3 emission had significantly different spectral properties when observed at Lick: Mrk 474, Mrk 728, and Mrk 1494 all appeared 
to be excellent reverberation mapping targets based on SDSS spectra observed in 2003-2006, but by 2011 the previously strong 
broa d H/3 lines faded so complete l y that they would be classified as Type 1.9 Seyferts. Transition s like these are not unprecedented 

(e.gJTohline & Osterbrockl 1 976llAretxaga et alJl999l:IShappee et alJ20~i4l:lDennev et alJ2014tlScott et al.ll2014HLaMassa et all 

120141) . but little statistical information is available on the timescale or duty cycle of such dramatic changes. For completeness, we 
present notes on the discarded targets below, along with plots of the spectra (Figure [24l) and comparisons with the earlier SDSS 
spectra for the three AGNs that exhibited a near-complete disappearance of broad H/3 (Figure [25l>. 

KA 1925+50: This Seyfert 1 galaxy at z = 0.067 is one of the brightest AGNs in the field observed during the main Kepler 
mission, and initial Kepler light curves were presented bv lMushotzkv et all (1201 li) . We confirm its Seyfert 1 classification, but its 
broad H/f emission line has a low amplitude relative to the host-galaxy continuum, so it is not an ideal target for H :’> rever beratio n 
mappi ng although Ha emission is very strong. Its full designation is 2MASX J19250215+5043134. See [Edelson & Malkanl 
(l2012h for further information about this object. 

Mrk 236: This AGN appeared to be in a fainter state than in an earlier SDSS spectrum, and was discarded in favor of other 
objects in the same RA range. 

Mrk 474 (NGC 5683): In the SDSS spectrum (observed 2003 April 4), this is a Seyfert 1 galaxy with very strong broad H/3 
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Fig. 25.— Comparison between SDSS and Kast spectra for Mrk 474, Mrk 728, and Mrk 1494, illustrating dramatic fading of the broad H/3 emission lines over 
a span of 5 to 8 years. The Lick spectra were not taken under photometric conditions, and have been scaled by an arbitrary factor to facilitate comparison. 


emission. Our Lick spectrum shows that the AGN has faded dramatically, and broad H/3 is no longer visible. 

Mrk 506: The broad H/3 emission in this galaxy is too weak to be useful for reverberation mapping, and our spectrum appears 
very similar to the SDSS spectrum. 

Mrk 728: This is another example of an object having strong broad H/3 in the SDSS spectrum (observed 2004 April 24) but 
showing little or no detectable broad H/3 in our Lick data. 

Mrk 885: This object is classified in NED as a Seyfert 1.0, but has relatively weak broad H/3 and is too faint and starlight- 
dominated to be a good reverberation-mapping target. 

Mrk 1494: The broad H/3 line in this AGN also appears to have faded dramatically in comparison with the SDSS spectrum, 
observed on 2006 May 28. The narrow [O III] and H/3 lines are also significantly fainter in the Lick spectrum than in the earlier 
SDSS data. 

NGC 6423 (Ark 524): This galaxy is listed in the lVeron-Cettv & Veronl d2006l) catalog and in NED as a Type 1.5 Seyfert, but 
we do not detect emission lines in its spectrum. 

B. ON THE RELATIONSHIP BETWEEN MEAN AND RMS PROFILE WIDTHS 

Broad-line widths measured from rms spectra are commonly used to derive black hole masses from reverberation data. In 
light of the very widespread use of rms profile w idths for determination of black hole masses and for calibrating the overall 
normalization of the AGN black hole mass scal e (lOn ken et al.1l2004t I Peterson et al.ll2004l : ICollin et al.ll2Q06b IWoo et al.ll20ldl : 

iGraham et alEoH iPark et al.ll2012bt;iGrier et al.ll2013hllHo&Kimll2014l) . we include in the following sections a description of 

simulations which may be helpful in clarifying certain aspects of the behavior of the rms profile in time-series spectroscopy of 
AGNs. 

We first examine the relationship between broad-line profile widths in mean and rms spectra. The mean line profile gives 
a useful representative view of the typical emission-line width during an observational campaign. However, the mean profile 
incorporates the complex and nonlinear response of the line profile t o con tinuum variations (including profile width variations 
and centroid shifts as well as flux variations, as illustrated in Section ROll . so it is not identical to the single-epoch profile that 
would be observed when the AGN is at its mean flux. The interpretation of the rms spectrum is less straightforward. It is 
sometimes stated that the rms s pectrum represents the variable portion of an emission line, because non varying components do 
not contribute to the rms (e.g.. IPeterson et afll 1 9981 fPeterson & WandeHl 19991 : iFromerth & MelialfeOOOl) . However, the BLR is 
not composed of distinct variable and nonvariable components; rather, different regions of the BLR have different responsivities 
to continuum variations, and the rms spectrum does not correspond directly to the spectrum emitted by any specific physical 
subcomponent of the BLR. Essentially, the rms spectrum is useful as a tool for visualization of the relative degree of variability as 
a function of wavelength across the line profile, due to spatial gradients in emission-line responsivity across the BLR combined 
with the effects of light-travel time across the BLR, but it is not obvious how to interpret the physical meaning of the velocity 
width of an rms line profile. Since reverberation-mapping campaigns are unable to acquire continuously sampled observations as 
an AGN varies, the detailed shape of both the mean and rms spectra depends additionally on the observational sampling pattern 
used during the campaign. Finally, continuum variability and noise in the data will have an impact on the line profiles and widths 
measured from rms spectra; we examine these effects in Appendix[C]below. _ 

The width of broad H/3 is typically f ound to be slightly narrower in rms profiles than in mean profiles (e.g. JSergeev et alll 1 999t 
ICollin et~aLll20Q6h iBentz et al.ll2010al) . Combining our sample with the previous 2008 Lick AGN Monitoring Project sample 
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(using measurements from lParket al.ll2012aL ba sed on spectral fitting method s similar to those employed here) plus additional 
Lick observations of Zw 229-015 (iBarth et al.l201 lal) and KA1858+4850 (lPei et al.l20l4i) . we have 17 independent measurements 
of the ratio cr(rms)/cr(mean) for H/3 in 15 AGNs, including two measurements each for Mrk 40 and Zw 229-015. This ratio ranges 
from 0.43 to 1.06, with a mean value of 0.85 and a_standard deviation of 0.15 among these 17 measurements. Similarly, in the 
sample of 15 AGNs examined bv lSergeev et al.l ( 119991) . the average value of the ratio cr(rms)/cr(mean) for the broad H/j line is 
0.87. 

The narrower widths of rms profiles for H/3 can be understood as a consequence of high-ionization gas in the inner BLR haying 
a lower responsivity to continuum variations than low-ionization gas more distant from the black hole; iKorista & Goadl ( 2004 ) 
predict a factor of ~ 5 change in H/3 responsivity across the radial range of the BLR. The high-velocity wings of the line profile 
are emitted primarily by low-responsivity portions of the BLR, while higher-responsivity gas contributes mainly to the line core. 
As a result, the core of the H/3 profile varies more strongly than the wings, making the rms profile narrower than the mean profile. 

Why is it that er(rms)/<T(mean) typically takes the value ~ 0.85? The purpose of this Appendix is to illustrate by means of 
simple numerical simulations that this ratio is dictated by the relationship between the broad H/3 single-epoch luminosity and 
line width (as described in Section^, which is fundamentally a consequence of the variation in emission-line responsivity across 
the BLR. The calculations presented here are not derived from a physical model of BLR variations; this is merely a numerical 
experiment designed to examine the behavior of rms profiles when the line width responds in a well-defined way to fl ux variatio ns. 

We begin with a Gaussian line profile as the baseline profile for steady-state emission. Using the method of iTimmer & Konigl 
(IT995h . we generate a random light curve of some specified duration, sampling, power spectrum shape, and rms variability 
amplitude (or F v ar ). As we describe below, the details of the light curve construction have little or no effect on the ratio of rms 
profile width to mean profile width. 

Then, we generate a time series of line profiles in which both the line luminosity and its width vary in response to continuum 
fluctuations. Assuming a delta function with lag r for the velocity-integrated transfer function, at each time t the integrated 
luminosity L of the emission line is made to scale linearly with the input continuum light curve luminosity at time t — T. The 
width of the Gaussian profile is also made to vary, such that the product cr 2 L 0 - 5 remains constant as the line luminosity varies. 
This represents the most common behavior seen in our sample [i.e., <r(H/3)oc /(H/3) -0 25 ], and corresponds to a situation in which 
the single-epoch virial product remains roughly constant as the AGN luminosity varies. From the time series of resultant line 
profiles, we construct the mean and rms profiles and calculate au ne for each. 

The result of this experiment is surprisingly simple. In the mean spectrum, the emission line has a profile that is very nearly 
Gaussian. As expected for a broad line showing higher variability amplitude in the core than in the wings, the rms profile is 
narrower than the mean profile. The rms profile has an approximately Gaussian shape, but exhibits additional structure in the 
faint high-velocity wings, and this profile shape is nearly independent of the input light curve properties. The peak height of the 
rms profile depends on the variability amplitude of the light curve, while its width is only very weakly dependent on variability 
amplitude. In the limit of small but nonzero light curve variance, the ratio (7(rms)/cr(mean) converges to 0.824. As the light 
curve variance is increased, we find that this ratio becomes very slightly sensitive to the input light curve properties. For example, 
when the input light curve is scaled to have F v ar = 0.2, the simulations produce a result of <r(rms)/a(mean) = 0.83 ±0.01, where 
the uncertainty represents the standard deviation of values from 1000 simulations carried out with different random light curves. 
Figure [26] shows an example of a mean and rms profile obtained from one such simulation with F mr = 0.2 over a light curve 
having 500 equally spaced points, with a fluctuation power spectrum of the form P{v) oc v~ 21 (where v is the temporal frequency 
of fluctuations). The rms profile shape would appear virtually identical for different choices of light curve duration or sampling, 
power spectrum slope, or variability amplitude, although its absolute normalization depends on f/ ar . 

These simulations produce a ratio cr(rms)/cr(mean) that is remarkably close to the observed mean ratio of 0.85 ±0.15 found 
from our Lick Observatory monitoring campaigns. In other words, the typical ratio of line widths observed in rms and mean line 
profiles simply reflects the fact that in most AGNs, the broad H/3 line maintains a roughly constant value of the product L il5 a 2 as 
the AGN luminosity varies. The larger dispersion of values observed in real AGNs most likely results from the range of different 
(non-Gaussian) profile shapes in different AGNs, and also differences in the scaling relationship between single-epoch luminosity 
and profile width. Changing the power-law relationship between line width and luminosity does alter the width of the rms profile 
relative to the mean profile. For example, if we calculate models for a flatter single-epoch scaling relationship of crii ne oc L 0,(17 
(as observed in Mrk 50), the simulations yield cr(rms)/er(mean) = 0.92 ±0.01. 

The prevailing use of rms profile widths as the basis for calibrating the AGN black hole mass scale underscores the need 
for further work to understand the rms spectrum by means of of detailed simulations of BLR line profiles incorporating BLR 
dyna mics a nd photoioni zati on physics, including luminosity-dependent changes in emission-line responsivity (Korista & Goad 
2004; Goad & Koristal2014l). The dynami cal modeling approach for interpretation of reverberation-mapping data dPancoast et ah 
2012 ; Li et al.l2013UPancoast et al.l2014al) avoids the problem of relying on simplistic measures of line width from either mean or 
rms spectra, and can provide more direct constraints on black hole masses. For future development of these models, it should be 
a high priority to incorporate a spatially variable and luminosity-dependent emission-line responsivity based on photoionization 
modeling of BLR structure, in order to more accurately model an observed time series of variable line profiles. If the responsivity 
is assumed to be spatially and/or temporally constant across the BLR, it will not generally be possible for these models to replicate 
the observed anticorrelation between line width and luminosity, or the relationship between mean and rms profile widths. 

C. BIASES IN RMS PROFILE WIDTHS 

Since line widths measured from rms spectra are used to derive black hole masses, it is useful to examine possible sources 
of systematic error in these measurements. The rms spectrum blends information about emission-line variations, continuum 
variations, and random noise in the data, and it can deviate substantially from the “ideal” rms spectrum that would be derived if a 
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Fig. 26.— Example from a simulation of a mean and rms line profile, constructed from a time series of Gaussian profiles having variable luminosities and 
widths as described in the text, with a 2 L 0,5 = constant. The amplitude of the rms profile is set by the overall variability amplitude of the input light curve in the 
simulation, but its shape is nearly independent of the light curve properties. 


single broad emission line could be observed without contamination from other spectral components or noise. In this Appendix, 
we describe simulations designed to examine the impact of continuum variability and photon-counting noise on line widths 
measured from rms spectra. 


Cl. Bias due to continuum variability 

This same simulation methodology described above allows us to examine a potential source of bias in the rms profile width that 
can occur when the rms spectrum is constructed from a set of spectra that have not had the variable AGN continuum removed. 
Components of constant flux, such as narrow emission lines or starlight, do not contribute to the rms spectrum except to the extent 
that they show some night-to-night variations due to seeing differences, miscentering of the AGN in the slit, random wavelength 
or flux calibration errors, or Poisson noise. However, it is not widely appreciated that the variable AGN continuum underlying 
an emission line can affect the shape of the rms profile of the line, and that this can potentially introduce significant biases into 
black hole virial masses determined from rms profile widths. 

The reason this bias occurs is that the continuum and emission-line variations are out of phase with one another: there will be 
portions of the light curves where the continuum is rising while the emission-line flux is falling, and vice versa. This out-of-phase 
variability alters the rms amplitude of the total (line plus continuum) flux, and tends to preferentially suppress the rms flux in the 
line profile wings relative to the line core. Consequently, the rms profile constructed from total-flux spectra can sometimes have 
a significantly different width than it would if the emission line could be observed in isolation without continuum contamination. 
Even a modest bias in determining the rms profile width can have a substantial and deleterious effect on the derivation of black 
hole masses, because black hole mass scales as the square of the line width. The typically employed procedure of fitting and 
subtracting the continuum from the total-flux rms spectrum prior to measuring the profile width does not alleviate this problem. 

To quantify this bias, we carried out simulations using Gaussian line profiles as described above. As before, the total line 
luminosity responds linearly to continuum variations, and the line luminosity and width vary such that a 2 L 0 5 remains constant. 
For the input light curves, we assumed a flu ctuatio n power spectru m of the form P(y) oc v~ 1J , similar to the properties of nearby 
AGNs monitored by the Kepler mission (iMushotzkv et al.ll201 ll) . and we carried out simulations with light curves of varying 
duration and various values of F var . The velocity-integrated transfer function was assumed to be a delta function at a lag of 10 
days. The emission line was set to have a mean equivalent width of 120 A relative to the mean AGN continuum flux density 
(typical for the broad H/3 lines in our sample) and an initial width of <j\i ne = 2000 km s -1 , and the continuum flux level was 
assumed to be flat over the wavelength range of interest. Spectra were constructed for velocity bins of 60 km s -1 , corresponding 
closely to the 1 A wavelength bins of the Kast blue-side data. These simulations do not include the effects of photon-counting 
noise, random wavelength or flux calibration errors, contributions of narrow-line or host-galaxy components to the spectra, or 
nonuniform sampling cadence, all of which can further alter the rms spectrum. (The impact of photon-counting noise is described 
in the next section.) 

Starting from each simulated continuum light curve, we generated a time series of emission-line profiles. At each time step, we 
created a spectrum of the emission line alone, and a total-flux spectrum including both the line and the continuum. Two versions 
of the rms spectrum were then computed from the time series of line-only and total-flux spectra. Then, the line dispersion 
fjline was measured from each of these. For the total-flux rms spectrum, the surrounding continuum level was removed prior to 
measuring the line dispersion by subtracting the minimum value of the rms spectrum. In order to minimize the impact of rms 
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Fig. 27. — Results from simulations to examine the ratio of line profile widths in total-flux rms spectra vs. emission-line only rms spectra. Left: Example of a 
simulation exhibiting a significant bias in the rms profile width when using the total-flux rms spectrum. The upper panels show the continuum and emission-line 
light curves; each point in the emission-line light curve corresponds to a simulated observation of the spectrum. The continuum has F var = 0.1, and the emission 
line has a lag of 10 days and an equivalent width of 120 A. The lower panel shows different versions of the rms spectrum. In black is the “true” rms spectrum 
determined from spectra of the emission line only (without continuum) over the 50 epochs of the simulation. The red dot-dashed curve shows the total-flux 
rms spectrum which includes the contribution of the variable continuum. The blue dotted curve shows the same total-flux rms spectrum, after subtracting the 
minimum continuum value and rescaling to match the peak flux of the true rms spectrum. In this case, the dispersion crij ne of the total-flux rms spectrum is 75% 
of the dispersion of the true rms spectrum, which would lead to an inferred black hole mass that was only 56% of its “true” value. Right: Simulation results 
showing the distribution of the ratio of emission-line dispersions crii ne measured from the total-flux and line-only spectra. Simulations were carried out for a lag 
of 10 days, monitoring durations of 30, 50, and 100 days with daily sampling, and F va r=0.1. 


flux in the very distant high-velocity line wings on the trii ne measurements, the dispersion was measured only over the wavelength 
range between the two minima on either side of the profile, in cases where the rms profile exhibited local minima straddling the 
line profile. Similarly, for the line-only rms spectrum we also measured the dispersion only over the wavelength range between 
the two local minima that occur in the profile wings. In real data, these very shallow minima would most likely not be visible 
as a result of the additional contributions to the continuum level of the rms spectra from other noise sources, and the rms profile 
width would be measured over some visually defined wavelength range corresponding roughly to the full width at zero intensity 
of the rms profile. 

The simulations show that the width of the total-flux rms line profile can often deviate very substantially from the width of 
the line-only rms profile, particularly when the monitoring duration is only a few times longer than the emission-line lag. Figure 
[27] (left panel) shows an example of one simulation, with a 50-day spectroscopic monitoring duration (i.e., 5 times longer than 
the lag, and assuming nightly observations with no gaps). The figure illustrates both the line-only rms profile, which should 
be considered the “true” rms profile for this simulation, and the total-flux (line plus continuum) rms profile. By subtracting the 
continuum level from the total-flux rms spectrum and rescaling it to match the peak of the true rms profile, it is immediately 
apparent that the total-flux rms profile is much narrower than the true profile. In this instance, er(total)/<T(true) = 0.75, which 
would bias the inferred black hole mass by a factor of 0.75 2 = 0.56. 

To examine the distribution of this bias factor, we carried out a set of these simulations for campaigns of duration 30, 50, and 
100 days, assuming nightly sampling in each case and a fixed value of F var = 0.1 for the AGN continuum. In Figure [27] (right 
panel), we plot the distribution of the ratio er(total)/cr(true) determined from a set of 10 4 simulations at each campaign duration. 
The bias distribution is a strong function of campaign duration. For example, we find that the percentage of simulations returning 
<r(total)/cr(true) < 0.75 is 22%, 6%, and 0.3% for campaign durations of 30, 50, and 100 days, respectively. In most of the recent 
large reverberation-mapping programs carried out by our team and others, the monitoring duration is usually at least ~ 5 times 
longer than the measured H/3 lags of the targets, but as future programs push toward higher-luminosity AGNs at higher redshifts 
there may be increasing numbers of AGNs for which the spectroscopic time series is sampled less intensively, over a monitoring 
duration that might be just barely sufficient to measure the lag (typically about ~ 3 times the lag). In that situation, the total-flux 
rms spectrum can frequently give a misleading measure of the emission-line width. 

In actual reverberation data sets, the magnitude of this bias will depend on the details of the light curve shape and variability 
amplitude, the transfer function shape, the profile shape and equivalent width of the emission line, and the sampling cadence of 
the observations, so the outcome will be different for every time series of data for every AGN. Our intention is not to carry out 
exhaustive simulations for a broad range of AGN parameters, but simply to demonstrate that this bias can occur under typical 
circumstances of a reverberation measurement when the total-flux rms spectrum is used. The simulation results shown in Figure 
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[27]illustrate that <7n ne can sometimes be severely underestimated even when the light curves have sufficient duration and sampling 
that the emission-line lag could easily be measured accurately. 

We propose that the best way to avoid this problem is to carry out continuum fitting and subtraction prior to construction of 
the rms spectrum, so that line widths can be measured from an rms spectrum that is free of continuum variability contamination. 
Other adva ntages of carrying out spectral decompositions for measurement of H/3 rms profile widths have been discussed by 
iPark et aD ( 2012al~) : they focused on the importance of removing He II, Fe II, narrow-line components, and starlight prior to 
constructing rms spectra, but they did not remove the nonstellar featureless continuum. Figure[6]illustrates the construction of an 
rms spectrum for the broad H/3 line in isolation, with other emission-line and continuum components removed, and in this paper 
we have measured the profile widths from rms spectra constructed in this way. It is still necessary to carry out a local continuum 
subtraction by fitting a linear continuum model to the rms flux on either side of the H/3 line in order to remove the small residual 
rms flux level due to random noise in the data. For the case of Mrk 40 illustrated in Figure [6] the values of <rii ne for the total-flux 
rms spectrum and the line-only rms spectrum differ by just 2%. Based on our simulations, only a very small difference would be 
expected given that the campaign duration was ~ 20 times longer than the H/3 lag of ~ 4 days for this AGN dBentz et al.ll2009bl) . 

C2. Bias due to photon statistics 

The finite S/N of spectroscopic observations is the source of a different bias that can afflict profile widths in rms spectra. The 
pixel values in an rms spectrum are determined not just by variability in the AGN, but also by the noise level in the spectra. A 
nonvarying object would still be observed to have a nonzero rms spectrum due to photon-counting noise and detector noise. The 
contribution of these noise sources can affect the profile width in the rms spectrum simply because the flux density is largest in 
the core of the emission line, hence the noise level will be highest there as well. The larger contribution of noise in the line core 
will elevate the peak of the line profile in the rms spectrum relative to the wings, making the line slightly narrower than it would 
appear in the absence of noise. While this is not a large effect, its impact is amplified because virial masses scale as the square of 
the line width. 

To examine the impact of S/N on rms profile widths, we carried out further simulations using the same methodology described 
above. Continuum light curves and the resulting time series of line profiles were created using the same prescriptions and 
parameters described in the previous section. The continuum and emission-line spectra for each simulation epoch were added 
to give the total flux spectrum at each epoch. Then, Gaussian random noise was added to each spectrum in the time series. The 
Gaussian dispersion for this pixel-to-pixel noise level was set by fixing the continuum S/N to a given value, and scaling the noise 
1 /2 

level by Sj' , where S, is the flux density at pixel i in the spectrum. Each simulation was run for a total spectroscopic monitoring 
duration of 100 days with nightly sampling, an emission-line lag of 10 days, and a specified value of F ml for the continuum 
variability amplitude. 

From each simulation, we computed a “true” rms spectrum from the time series of noise-free emission-line spectra only, with 
no continuum contributions. Using the noise-added data, we subtracted the continuum flux level from each nightly spectrum 
and constructed a line-only rms spectrum that included the noise contribution. The line dispersion of the true rms profile was 
measured as described previously, by including the flux only between the two local minima on either side of the rms profile. To 
measure the width of the noise-added rms profile, we fit and subtracted the fiat noise continuum using a straight-line fit, and then 
calculated the line dispersion over the same velocity range that was used for the true rms profile measurement. The ratio of the 
noise-added rms line width to the true rms line width gives the bias factor, cr(noise-added)/(j(noise-free). Since the continuum 
flux was subtracted from each noise-added spectrum prior to constructing the noise-added rms spectrum, this measurement is 
immune to the continuum-variability bias described in the previous section. 

The degree of bias in the rms profile width depends primarily on F var and on the S/N of the data. We carried out simulations for 
continuum S/N ranging from 20 to 100 (in increments of 2) and for F var ranging from 0.05 to 0.30 (in increments of 0.01). For 
each set of these parameters, we carried out 1000 simulation trials and determined the median bias factor. Figure[28]displays the 
results, illustrating the ratio cr(noise-added) / at noise-free) as a function of F var and the continuum S/N. We also display the bias 
in the virial product, VP(noise-added)/VP(noise-free) =[cr(noise-added)/(r(noise-free)] 2 . For typical parameters of a successful 
reverberation measurement, corresponding to F var >0.1 and S/N of 30-100 per pixel, the rms profile width is biased low by 
~ 3-7%, and virial products are biased low by ~ 6 - 13%. The bias vanishes in the limit of high S/N and high variability 
amplitude, but even at S/N=100 per pixel the virial products are biased low by ~ 5 - 10%. Unlike the continuum-variability bias 
described in the previous section, this effect is not alleviated by having a long and well-sampled campaign with a large number of 
observations. For given values of S/N and F var , increasing the number of spectra in the time series will not remove this systematic 
effect of photon-counting noise on the rms profile width, although it will make the rms spectrum less noisy. 

Given the data quality for typical reverberation campaigns, these simulations suggest that virial products calculated using rms 
spectrum profile widths will be systematically low typically by ~ 10%. Th e sta ndard practice o f normalizing reverberation- 
based black hole masses to the local Mbh - cr relation of inactive galaxies (lOnken et al.ll2004t IWoo et al. 1201 Oh will compensate 
for this offset, and consequently the the virial /-factor derived from this procedure (e.g.. lOnken et all 20041 IWoo et al.ll2010t) 
will be biased high by ~ 10%. The impact of this effect is much smaller than the overall uncertainty on the mean / factor 
for the set of reverberation-mapped AGNs, so this does not currently represent a major problem for derivation of AGN black 
hole masses, but nevertheless it is a real effect that is present in all reverberation-mapping data. IPark et afl (12012a!) presented a 
maximum-likelihood method to calculate rms spectra that was designed to be immune to this problem; further tests of this or 
other alternative methods to calculate rms spectra would be worthwhile. 
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